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Abstract 



Quantum chaos is the study of quantum systems whose classical description is 
chaotic. How does chaos manifest itself in the quantum world? In this spirit, we 
study the dynamical generation of entanglement as a signature of chaos in a system 
of periodically kicked coupled-tops, where chaos and entanglement arise from the 
same physical mechanism. The long-time entanglement as a function of the position 
of an initially localized wave packet very closely correlates with the classical phase 
space surface of section - it is nearly uniform in the chaotic sea, and reproduces 
the detailed structure of the regular islands. The uniform value in the chaotic sea 
is explained by the random state conjecture. As classically chaotic dynamics take 
localized distributions in phase space to random distributions, quantized versions 
take localized coherent states to pseudo-random states in Hilbert space. Such ran- 
dom states are highly entangled, with an average value near that of the maximally 
entangled state. For a map with global chaos, we derive that value based on new an- 
alytic results for the entropy of random states. For a mixed phase space, we use the 
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Percival conjecture to identify a "chaotic subspace" of the Hilbert space. The typical 
entanglement, averaged over the unitarily invariant Haar measure in this subspace, 
agrees with the long-time averaged entanglement for initial states in the chaotic sea. 
In all cases the dynamically generated entanglement is that of a random complex 
vector, even though the system is time- reversal invariant, and the Floquet operator 
is a member of the circular orthogonal ensemble. 

Continuing on our journey to find the footprints of chaos in the quantum world, 
we explore quantum signatures of classical chaos by studying the rate of information 
gain in quantum tomography. The measurement record is obtained as a sequence 
of expectation values of a Hermitian operator evolving under repeated application 
of the Floquet operator of the quantum kicked top on a large ensemble of identical 
systems. We find an increase in the rate of information gain and hence higher 
fidelities in the process when the Floquet maps employed increase in chaoticity. We 
make predictions for the information gain using random matrix theory in the fully 
chaotic regime and show a remarkable agreement between the two. Finally we discuss 
how this approach can be used in general as a benchmark for information gain in an 
experimental implementation based on nonlinear dynamics of atomic spins measured 
weakly by the Faraday rotation of a laser probe. 

The last part of this thesis is devoted to the study of the nature of quantum 
correlations themselves. Quantum correlations are at the heart of the weirdness of 
quantum mechanics and at the same time serve as a resource for the potential bene- 
fits quantum information processing might provide. For example, Einstein described 
quantum entanglement as "spooky action at a distance" [1]. However, even entangle- 
ment does not fully capture the complete quantum character of a system. Quantum 
discord aims to fill this gap and captures essentially all the quantum correlations in a 
quantum state [2] . There is a considerable interest in the research community about 
quantum discord, since there is evidence showing this very quantity as responsible 
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for the exponential speed up of a certain class of quantum algorithms over classical 
ones [3]. 

Now, an important question arises: Is discord just a mathematical construct or 
does it have a definable physical role in information processing? This thesis provides 
a link between quantum discord and an actual physical task involving communication 
between two parties [ 1] . We present an operational interpretation of quantum discord 
based on the quantum state merging protocol. Quantum discord is the markup in the 
cost of quantum communication in the process of quantum state merging, if one dis- 
cards relevant prior information. We further derive a quantitative relation between 
the yield of the fully quantum Slepian-Wolf protocol in the presence of noise and the 
quantum discord of the state involved. This protocol is the most general known in 
the family of protocols in quantum information theory, a unification of essentially 
all bipartite, unidirectional and memoryless quantum communication protocols. The 
significance of quantum discord in noisy versions of teleportation, super-dense cod- 
ing, entanglement distillation and quantum state merging are discussed. We also 
demonstrate similar roles for quantum discord in quantum computation and corre- 
lation erasure. Our work shows that quantum discord captures and quantifies the 
advantage of quantum coherence in quantum communication. 
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Chapter 1 
Introduction 



Information is physical. As first emphasized by Landauer [5], this indicates that 
the information processing capabilities of a device are not separate from the physics 
that governs its operation. On the one hand, over the last few decades, the field 
of quantum information science has shown that the devices governed by the laws of 
quantum mechanics can have information processing capabilities superior to their 
classical counterparts. It has also shed light on what properties of quantum systems 
should be harnessed for building future technology and engineering applications. 
Here, the role of information, quantified suitably as entanglement and other measures 
of quantum correlations between sub-systems, provides vital clues to the superior 
information processing capabilities of devices based on quantum mechanics. 

On the other hand, and at a more fundamental level, understanding how physical 
systems process and exchange information is crucial to gaining insights into the work- 
ings of our universe. For example, the connections between entropy, information and 
thermodynamics form the cornerstone of statistical mechanics. An example more 
relevant to us is Feynman's path-integral formulation of non-relativistic quantum 
mechanics [6]. In this approach, physical phenomena are described by events. An 
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event can occur through various alternatives/paths, each of which is characterized 
by a complex probability amplitude, that has both real and imaginary components 
in general. If these paths are in principle indistinguishable, i.e. there is no informa- 
tion whatsoever in the universe that can help us distinguish between them, then in 
that case the corresponding probability amplitudes add up causing interference. For 
example, consider the 2-slit experiment. As long as there is no information available 
as to which slit the photon takes, we see interference fringes on the screen. The 
photon is considered to be in a superposition of two wave packets, centered around 
the classical path out of each slit. The question as to why does the universe care 
about which "path" the photon takes is both philosophical and intriguing. It cer- 
tainly tells us that the universe cares about certain kinds of information or the lack 
of it in formulating its laws. 

We saw above that a photon is considered to be in a superposition of two wave 
packets, centered around the classical path out of each slit. In general, the system is 
described by a vector in the Hilbert space. This vector in general is a superposition 
of N basis vectors (this comes from the linearity of Schrodinger's equation), where 
TV is the dimension of the Hilbert space. For an infinite dimensional Hilbert space, 
the sum in the superposition is replaced by an integral. The macroscopic world 
around us does not permit the majority of these state vectors. For example, a tennis 
ball is never seen in a superposition of two values of position simultaneously. This 
leads us to another fundamental question : How does the classical physics governed 
by Newton's laws emerge out of the underlying quantum mechanics? The flow of 
information in quantum physics is at the heart of understanding the emergence of 
the classical world and the disappearance of quantum coherence via decoherence. 
Decoherence, which is the leaking of information to the environment, suppresses 
quantum interference effects and enables the quantum to classical transition. Thus 
we see that information plays a fundamental role, perhaps as fundamental as that of 
energy and momentum, in the workings of our universe. 
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Chapter 1. Introduction 

The connections between information and dynamical properties of quantum sys- 
tems is a recurring theme of this thesis. The first half of this thesis looks at quantum 
chaos from a quantum information science perspective. 

1.1 Classical and quantum chaos 

Broadly speaking, classical chaos is the long term aperiodic behavior in a determin- 
istic dynamical system that exhibits sensitive dependence on initial conditions [7]. 
Aperiodic long term behavior means that there are trajectories which do not settle 
down to fixed points, periodic orbits, or quasiperiodic orbits in the limit t — > oo, 
where t is the time of evolution of the trajectory. Sensitive dependence on initial 
conditions means that nearby trajectories separate exponentially fast; the rate of 
separation is given by the Lyapunov exponent, A, which characterizes the dynamics 
of the system. A system with N degrees of freedom having N constants of motion 
is said to be integrable and its dynamics is regular. When one has fewer than N 
constants of motion, then the individual trajectories can explore the phase space in 
a complex manner and the system can exhibit chaos. 

It is not difficult to see that the above "definition" of chaos fails in the quantum 
domain. A quantum state is not a point in the phase space but is described by a 
state vector. The time evolution of the state vector due to the Schrodinger's equation 
is unitary. This means that the overlap of two state vectors undergoing evolution 
is constant with time. Furthermore, while classical chaos can lead to infinitely fine 
structures in the phase space, in quantum mechanics, Planck's constant, h, sets the 
scale for such structures. According to Heisenberg's uncertainty principle, our res- 
olution of the phase space is determined by Planck's constant. This is often stated 
as the key reason for the absence of chaos in the quantum domain. This however is 
not the complete story. An alternate description of classical mechanics, involving the 
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evolution of classical probability densities, preserves the distance between two prob- 
ability densities as a function of time. Hence, the distance between two probability 
densities does not show exponential sensitivity even for classical mechanics. 

All this leads to two interesting questions: 

1. How does classically chaotic dynamics inform us about certain properties of 
quantum systems, e.g. the energy spectrum, nature of eigenstates, correlation 
functions, and more recently, entanglement. Alternatively, what features of 
quantum systems arise due to the fact that their classical description is chaotic? 

2. Since all systems are fundamentally quantum, how does the classical chaos, 
with trajectories sensitive to initial conditions, arise out of the underlying quan- 
tum equations of motion? 

These two questions are not unrelated. However, the first question deals mainly with 
finding the signatures of chaos by studying the properties of the quantum Hamilto- 
nian, while the second concerns with the dynamical behaviour of quantum states and 
the emergence of classically chaotic behaviour. 

One of the first steps in addressing question (1) was taken by Gutzwiller [8] by 
developing the trace formula. The "trace formula" relates the density of states of 
the quantum system in the semi-classical regime to the properties of the classical 
periodic orbits. A central result of quantum chaos is the connection with the the- 
ory of random matrices [9]. In the limit of large Hilbert space dimensions (small 
h), for parameters such that the classical description of the dynamics shows global 
chaos, the eigenstates and eigenvalues of the quantum dynamics have the statistical 
properties of an ensemble of random matrices. The appropriate ensemble depends 
on the properties of the quantum system under time-reversal [9]. In this direction 
an important quantum signature of chaos was obtained by Bohigas and collabora- 
tors [i l], describing the spectral statistics of quantum Hamiltonians whose classical 
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counterparts exhibit complete chaos using random matrix theory. Such signatures of 
quantum chaos have mainly focussed on the time-independent Schrodinger's equation 
and features like energy spectra and eigenstates. The connections between quantum 
chaos, time reversal symmetry and Random Matrix Theory will be important to us, 
and hence we devote the next chapter giving an overview of these aspects. 

In studying the dynamical aspects of quantum chaos, it is worth mentioning the 
conflict between the classical and quantum theories in describing the future dynamics 
of system. Ehrenfest's correspondence principle gives us an estimate of how well the 
quantum expectation values follow Hamilton's canonical equations of motion. The 
"quantum break time" is the characteristic time at which the two diverge. This 
correspondence follows when the width of a well localized quantum state is extremely 
narrow as compared to the characteristic variations in the macroscopic potential. 
For systems having regular dynamics, significant deviations between the quantum 
expectation values and classical dynamics occur on very late time scales [13, 14]. 
These deviations scale roughly with the size of the system in action relative to h. 
This, however, is not true for chaotic systems. In the case of chaotic systems, the 
break time can be logarithmic in the size of the action relative to h, and thus even for 
macroscopic systems, this time will be surprisingly short. Zurek took this argument 
to its absurd conclusion and showed how the trajectory of Hyperion, a moon of 
Saturn, will diverge from classical behavior after only 20 years! He concluded that 
decoherence resolves the paradox [2, 16]. 

Habib and collaborators picked up this problem almost a decade ago and stud- 
ied it from the point of view of continuous measurement and quantum trajectories 
[12]. In particular, they asked the following question: Under what conditions does a 
quantum trajectory that tracks a measurement of a given observable follow the clas- 
sical trajectory? For the case of a single degree of freedom they showed conditions 
under which this was true. In particular they showed that if the backaction was 
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strong enough that the wave packet remain sufficiently localized that interference 
across different classical trajectories is negligible, but weak enough that the backac- 
tion noise is small compared to the quantum uncertainty, then the system follows 
essentially classical trajectories. These conditions could be satisfied if the system 
is sufficiently macroscopic and the measurement backaction is sufficiently strong. 
Ghose and collaborators extended this to study the quantum to classical transition 
in a bipartite system consisting of coupling between spin and motional degrees of 
freedom [17]. They concluded that even when the measured motional degree of free- 
dom can be treated classically, the entanglement between the two subsystems causes 
strong measurement backaction on the quantum spin subsystem and one does not 
recover classical trajectories in this regime. It is important for both subsystems to 
be sufficiently macroscopic for the recovery of classical trajectories. Habib et al. 
extended this to study the Lyapunov exponents of quantum trajectories even when 
the system was not macroscopic [18]. 

Though quantum systems show no exponential separation under the evolution 
of a known unitary evolution, they do show a sensitivity to the parameters in the 
Hamiltonian [] ]. Peres [19] showed that the evolution of a quantum state is altered 
when a small perturbation is added to the Hamiltonian. As time progresses, the 
overlap of the perturbed and unperturbed states gives an indication of the stability 
of quantum motion. It was shown that if a quantum system has a classically chaotic 
analog, this overlap has a very small value. On the other hand, if the classical analog 
is regular, the overlap remains appreciable. In another perspective, as seen in the 
work of Schack and Caves, quantum systems exhibit chaos when they are perturbed 
by the environment. They become hypersensitive to perturbations [ ], as seen in 
the information theoretic studies of the cost to maintain low entropy in the face of 
loss of information to the environment. This particular feature of quantum chaotic 
systems has several interesting consequences. For example, Shepelyansky has done 
extensive work on the issue of many-body quantum chaos in the quantum computer 
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hardware and its effect on the accuracy of quantum computation [23]. Recently 
classical simulations of quantum dynamics have been connected to integrability and 
chaos [25]. 

It is imperative to mention the role played by quantum information theory in 
the above journey. Quantum information science has added a whole new perspective 
to the study of quantum mechanics. This has resulted in a better understanding of 
quantum phenomena like entanglement and decoherence, and given us the tools to 
view certain quantum properties of physical systems as a resource. This has also 
enabled us to address the key questions in quantum chaos from a new perspective. 
As mentioned above, this has led to an information theoretic characterization of 
quantum chaos [20] and enabled the exploration of the behavior of chaotic quantum 
systems in the presence of environment induced decoherence [26] along with its con- 
nection to the quantum to classical transition. The study of quantum chaos from a 
quantum information perspective is also closely related to the theory and application 
of random quantum circuits [24]. 

A major part of this thesis is devoted to the study of signatures of quantum 
chaos as seen in the light of quantum information theory. We explore the connections 
between entanglement, information generation in quantum tomography, and chaos. 
Our studies throw light on the connections between ergodicity, entropy production 
and chaos. Moreover, employing continuous measurement quantum tomography, we 
develop an entirely novel paradigm to study quantum chaos. 

Here I introduce and motivate the work undertaken in this thesis in the above 
quest. I also mention the chapters where the reader can find the details of the 
projects undertaken. 
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1.2 Dynamically generated entanglement as a sig- 
nature of quantum chaos 

The connections between complexity, nonlinear dynamics, ergodicity, and entropy 
production, have long been at the heart of the foundations of statistical physics. A 
central goal of "quantum chaos" has been to extend this foundation to the quan- 
tum world. In the last decade, the tensor product structure of quantum mechanics, 
essential for understanding systems with multiple degrees of freedom, has come to 
the fore. In that context, one is naturally led to consider how the dynamical gen- 
eration of entanglement between quantum subsystems is connected with the chaotic 
dynamics of coupled classical degrees of freedom. Such studies address fundamental 
issues of complexity in quantum systems and are potentially applicable in quantum 
information processing, where entanglement is considered to be an essential resource. 

The connection between chaos in the classical description of Hamiltonian dynam- 
ics and entanglement in the quantum description has been the subject of extensive 
study over the last decade. The original motivation of Zurek and Paz was to ad- 
dress the quantum-to-classical transition [26]. By conjecturing that chaotic systems 
decohere exponentially fast through their entanglement with the environment, they 
hoped to resolve a paradox in which a macroscopic system would exhibit the effects 
of quantum coherence on a time scale logarithmic in h. 

To address these issues, we consider a model system of kicked coupled-tops, de- 
scribed in detail in Chapter 3. This system is motivated by its connection to possible 
experimental realizations, our ability to easily visualize the classical phases space, 
and to analyze the Floquet map. We use this system as a forum to explore ques- 
tions (1) and (2) above - how is chaos in the classical dynamics of the joint system 
correlated with the long-time averaged dynamically generated entanglement? 
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1.3 Quantum signatures of chaos in quantum to- 
mography 

At a fundamental level chaos represents "unpredictability", so this seems at odds 
with the goal of gaining information in order to estimate an unknown quantum state. 
However, on the flip side, this unpredictability represents the potential "information 
to be gained" in an estimation process. If everything is predicted and known, we learn 
nothing new. The missing information in deterministic chaos is the initial condition. 
A time history of a trajectory at discrete times is an archive of information about 
the initial conditions given perfect knowledge about the dynamics. Moreover, if 
the dynamics is chaotic the rate at which we learn information increases due the 
rapid Lyapunov divergence of distinguishable trajectories and we expect unbiased 
information because of the ergodic mixing of phase space. That is, if the information 
is generated by chaotic dynamics, the trajectory is random, and all initial conditions 
are equally likely until we invert the data and discover the initial state. 

The standard way to perform quantum tomography is to make projective mea- 
surements of an "informationally complete" set of observables and repeat them many 
times. The statistics obtained are used to estimate the expectation values of the ob- 
servables and hence the unknown initial state. The projective measurements pose 
a hurdle in exploring the connections between information gain in tomography and 
chaos due to large measurement back-action on the system. However, we overcome 
this by employing the protocol for tomography via weak continuous measurement 
developed by Silberfarb et al. 7]. In this protocol, the ensemble is collectively 
controlled and probed in a time dependent manner to obtain an "informationally 
complete" continuous measurement record. We consider the case of a very weak 
measurement such that the back-action is negligible. Our protocol for quantum to- 
mography via continuous measurement of a driven system [ '] gives us a window 
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into the complexity of quantum dynamics and its relationship to chaos. This work 
is intimately related to the protocols that have recently been implemented in the 
laboratory [28]. 

As in the case of deterministic chaos, the missing information in tomography is the 
initial condition. We try to accurately model all of the quantum dynamics occurring 
in the system, and then use the measurement time history to give us information 
about the initial quantum state. The dynamics is "informationally complete" if the 
time history contains information about an arbitrary initial condition. Our goal 
is to characterize and quantify the performance of tomography, when the dynamics 
driving the system are chaotic in the classical limit. We use this to draw comparisons 
between the role played by regular and chaotic dynamics in the information gain in 
the tomography procedure. In chapter 4, we report the novel and intriguing quantum 
signatures of classical chaos we have under this paradigm. 

The second half of the thesis is concerned with the study of quantum correlations 
and exploring their role in quantum information processing protocols. 

1.4 Nonclassical correlations and their role in 
quantum communication 

Quantum information science is primarily aimed at harnessing the quantum struc- 
ture of nature for information processing and computing tasks [29]. This quest has 
met with considerable success over the last 20 years, but there has been substantial 
progress in the other direction as well. Information theory has provided a novel 
framework for unraveling the intricacies of quantum mechanics. Quantum correla- 
tions, as well as classical ones are now viewed as resources, whose interconvertibility 
is governed by quantum information theory [30]. Foremost amongst these is evidently 
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entanglement, which provides enhanced performance in several important tasks like 
communication, computation, metrology and others [31]. However, even entangle- 
ment does not fully capture the complete quantum character of a system. 
Quantum discord aims to fill this gap and captures essentially all the quantum cor- 
relations in a quantum state [2]. There is a considerable interest in understanding 
the meaning of quantum discord, since there is evidence showing this very quantity 
as responsible for the exponential speed up of a certain class of quantum algorithms 
over classical ones [3]. Now, an important question arises: Is discord just a mathe- 
matical construct or does it have a definite physical role in information processing? 
This thesis provides a link between quantum discord and an actual physical task 
involving communication between two parties [I]. In chapter 5, we present an op- 
erational interpretation of quantum discord based on the quantum state merging 
protocol. Quantum discord is the markup in the cost of quantum communication in 
the process of quantum state merging, if one discards relevant prior information. 
In chapter 6, we further derive a quantitative relation between the yield of the fully 
quantum Slepian-Wolf protocol [33] in the presence of noise and the quantum dis- 
cord of the state involved. This protocol is the most general known in the family 
of protocols in quantum information theory, a unification of essentially all bipartite, 
unidirectional and memoryless quantum communication protocols. The significance 
of quantum discord in noisy versions of teleportation, super-dense coding, entangle- 
ment distillation and quantum state merging are discussed. We also demonstrate 
similar roles for quantum discord in quantum computation and correlation erasure. 
Our work shows that quantum discord captures and quantifies the advantage of 
quantum coherence in quantum communication. 
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1.5 Related Work and List of Publications 

Several new results are presented in this thesis. The work presented here has been 
published in refereed journals and conferences. 

1. Entanglement and generation of random states in the quantum chaotic dy- 
namics of kicked coupled tops, Collin M. Trail, Vaibhav Madhok, and Ivan H. 
Deutsch, Phys. Rev. E, 76, 046211, (2008). (Chapter 3) 

2. Quantum signatures of chaos in quantum tomography, Vaibhav Madhok, Carlos 
Riofrio, and Ivan H. Deutsch in preparation. (Chapter 4) 

3. Interpreting quantum discord through quantum state merging, Vaibhav Mad- 
hok and Animesh Datta, Phys. Rev. A 83, 032323, (2011). (Chapter 5) 

4. Role of quantum discord in quantum communication, Vaibhav Madhok and 
Animesh Datta, arxiv:1107.0994 (2011), 6th Conference on Theory of Quantum 
Computation, Communication and Cryptography ( Springer's Lecture Notes 
in Computer Science, 2011). (Chapter 6) 

5. Quantum discord as a resource in quantum protocols, Vaibhav Madhok and 
Animesh Datta, submitted to Phys. Rev. Lett. (Chapter 6) 

6. Quantum discord : The advantage of coherence in quantum communication, 
Vaibhav Madhok and Animesh Datta invited article to appear in the Interna- 
tional Journal of Modern Physics B, (Chapter 2 and Chapter 6). 

In addition, I have worked on quantum metrology during my graduate career at 
the UNM. The goal here is to study the fundamental bounds imposed upon precision 
measurements by quantum mechanics and to investigate the nature of quantum states 
and resources that might help us attain those bounds. In particular, we showed that 
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mode entangled, two mode, Gaussian states can be used to estimate a phase with 
a measurement uncertainty 5<p that scales as 1/N, where N is the average photon 
number in the state. We also considered the effect of photon loss on the measurement 
uncertainty. The results of this investigation will be reported elsewhere. 
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The purpose of this chapter is to build a foundation to aid understanding of the 
original work presented in the remaining part of this thesis. This chapter gives 
a pedagogical treatment of some technical aspects of this research, which will be 
helpful for understanding the work done and the results in the subsequent chapters. 

Broadly, studies in quantum chaos are divided into kinematics and dynamics. 
The cornerstone of the former, denoted by Berry as "quantum chaology" , is random 
matrix theory [ ], connecting the statistics of the spectrum and eigenvectors of 
Hamiltonians and Floquet maps of chaotic systems with those of random matrix 
ensembles. Classic works on the subject including level statistics [46], properties of 
Wigner functions [47], and quantum scars in ergodic phase spaces [55] have tended 
to focus on the properties of wave mechanics, e.g. the dynamics of single particle 
billiards [48] (also seen in the properties of classical waves, e.g. microwave cavities 
[49]). 

The connections between quantum chaos, time reversal symmetry and the random 
matrix theory will be important to us, and hence we start by reviewing these aspects 
of our study. 
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2.1 Quantum chaos 
2.1.1 Time reversal 

Roughly speaking, if one considers the time evolution of a physical system to be a 
movie being played, then the absence/presence of the time reversal symmetry is our 
ability /inability to tell which way, forward or backwards, the movie is running. For 
classical systems, time reversal symmetry implies that for any given solution x(£), 
p(t) of Hamilton's equations an independent solution x'(t'), p'(t') can be obtained 
with i! = —t and some algebraic operation relating x' and p' to the position x and 
momenta p [ ]. It is important to mention the distinction between irreversibility 
and the absence of time reversal invariance symmetry. Reversibility is related to 
whether or not we have lost information, and whether the past condition can be 
recovered from the present. One can have reversibility without having time reversal 
invariance symmetry. For example, the motion of a charged particle in a magnetic 
field is not time-reversal invariant but it is reversible since there is no dissipation or 
loss of information. 

In quantum mechanics, the quantum system obeying Schrodinger's equation 



is time reversal invariant, if for any given solution ip('x,t), there is another solution, 
,t'), with t' = —t and ip' uniquely related to ip [')]. For example, consider a 
spinless particle with the Hamiltonian described solely by a scalar potential 



If ^(x, t) is a solution satisfying the Schrodinger's equation with time t then, for 
t' = —t, , 0*(x, t') is a solution. Hence i/j ; (x,t) = Kip{yL,—i) solves the Schrodinger's 
equation where K is the complex conjugation operator defined with respect to the 



(2.1) 




(2.2) 
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position basis. In general, a complex conjugation operator K' can be defined with 
respect to any representation. 

Wigner showed that all time reversal operators T must be antiunitary [ 

<7W> = (2.3) 

Therefore, any time reversal operator T can be given the standard form 

T = UK (2.4) 

where U is a suitable unitary operator and K is the complex conjugation with respect 
to a standard representation. Another requirement for T is that any quantum state 
should be reproduced, to within a phase factor, on application of T twice, i.e 

T 2 = a, \a\ = l. (2.5) 
Using Eq. 2.4 and Eq. 2.5 one gets a useful condition on T, 

T 2 = ±1. (2.6) 
When we have a time reversal symmetry, and the antiunitary operator T is such that 

[H,T] = 0,T 2 = 1, (2.7) 

then the Hamiltonian can always be given a real matrix representation and such a rep- 
resentation can be found without diagonalizing H. In such a scenario, a T-invariant 
basis, {ipi,ip2, ■■■ipn} can be constructed such that Tipi = ipx. With respect so this 
basis, the Hamiltonian H = THT is a real matrix [9]. The canonical transformations 
for such Hamiltonians are the orthogonal matrices O, 00 T = 1. 

In the field of quantum chaos, time reversal becomes important in the following 
way. In the limit of large Hilbert space dimensions (small h), for parameters such 
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that the classical description of the dynamics shows global chaos, the eigenstates and 
eigenvalues of the quantum dynamics have the statistical properties of an ensemble 
of random matrices. The appropriate ensemble depends on the properties of the 
quantum system under time- reversal [9]. The ensemble of random matrices used to 
describe the Hamiltonians unrestricted by the time reversal symmetry is the Gaus- 
sian Unitary Ensemble (GUE). Similarly, the ensemble of random matrices used to 
describe the Hamiltonians having a time reversal symmetry are given by the Gaus- 
sian Orthogonal Ensemble (GOE). Since the random matrix theory is crucial to our 
understanding of the quantum chaos, we take look at it in the next section. 



2.1.2 Random Matrix Theory 



In 1984, Bohigas, Giannoni, and Schmit (BGS) conjectured that random matrix the- 
ory (RMT), originally introduced by Wigner and Dyson to explain the spectra of the 
complex many-body physics of nuclei, applied also to few body (even one body with 
hard wall boundaries) systems when the dynamics is completely chaotic. Using RMT 
to describe the properties of eigenvalues/eigenstates in completely chaotic system is 
the foundation of much of modern "quantum chaology" . There is overwhelming ev- 
idence for the existence of universality in local fluctuations in the quantum energy 
spectra for systems that display global chaos in their classical phase spaces. All 
such Hamiltonian matrices of sufficiently large dimension have spectral fluctuations 
which can be described by studying such fluctuations for an appropriate ensemble of 
random matrices [9]. Any member of such an ensemble can serve as a model of the 
Hamiltonian. 

In the field of quantum chaos, we typically look at twok kinds of random matrices 
- random Hermitian matrices and random unitary matrices. Random Hermitians 
are of interest when one wants to understand energy levels of a complex system. 
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The ensemble of random matrices used to describe Hamiltonians unrestricted by 
the time reversal symmetry is the Gaussian Unitary Ensemble (GUE). Similarly, 
the ensemble of random matrices used to describe those Hamiltonians having a time 
reversal symmetry are given by the Gaussian Orthogonal Ensemble (GOE). The other 
class of random matrices typically studied are the random unitary matrices. They 
are employed for periodically driven systems, as models of the unitary "Floquet" 
operators, F, describing the change of the quantum state during one cycle of the 
driving. Powers of the "Floquet" operator, F n , gives us a stroboscopic description 
of the dynamics. The ensemble of random unitaries are also known as the "circular 
ensembles" , originally introduced by Dysan. As was the case for random Hermitian 
matrices, time reversal symmetry arguments play a similar role in the choice of 
the appropriate ensemble of random unitaries employed to model the "Floquet" 
operator to study the properties of the chaotic system. Depending on whether the 
system has time reversal symmetry or not, the appropriate ensemble of random 
unitaries is called the Circular Orthogonal Ensemble (COE) or the Circular Unitary 
Ensemble (CUE) respectively. The eigenvectors of the COE and CUE have the 
same properties as that for the respective GOE and GUE, but the eigenvalues are 
distributed differently. The circular unitary ensemble (CUE) is just the ensemble 
of random unitary matrices picked from U(n) according to the Haar measure. CUE 
eigenvalues lie on the unit circle in the complex plane, hence the name. We look at 
an example of the construction of a Gaussian ensemble of Hermitian matrices below. 



An example of Gaussian Ensemble of Hermitian Matrices 

To illustrate the theory discussed above, we provide an example of construction of the 
Gaussian ensemble by constructing real symmetric 2x2 matrices with the orthogonal 
group as their group of canonical transformations as discussed in [9]. Consider the 
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matrix, 

Hu H12 
H 2 i H 22 

We want to determine the probability density, P{H), for the three independent 
matrix elements Hu, if 22, and H± 2 normalised as 

/+00 

P(H)dH n dH 22 dH 12 = 1 (2.8) 
-00 

To determine P(H) uniquely, we demand 

(1) P(H) must be invariant under any canonical transformation, i.e., 

P(H) = P(H'), H' = T HO, T = 0-\ (2.9) 

(2) The three independent elements must be uncorrelated. This implies that 
P(H) be of the form 



O 



P(H) = Pu(H n )P22(H 22 )P 12 (H 12 ). 
Considering infinitesimal orthogonal change of basis, 

1 9 
-6 1 

for which H' = T HO gives 

H'u = Hu — 26Hi 2 
H 22 = H 22 + 26H\ 2 
H[ 2 = H 12 + Q{H 12 — H 22 ). 



(2.10) 



(2.11a) 
(2.11b) 
(2.11c) 



Using the above equations, for infinitesimal 6, Pu(H' u ) = Pu(H n ) — 29H 12 dP d ^ L \ 
with similar expressions for P 22 (H' 22 ) and Pi 2 (H[ 2 ). Since P{H) is invariant under 
the orthogonal transformation we get, 

P(H) = P(H){1 - 9[2H 12 d -^ - 2H 12 ^ - (H U - H 22 ) d -^}}. (2.12) 



dH 



11 



dH 



22 



dH 



12 
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Since the angle 6 is arbitrary, its coefficient in the above equation must vanish, 

1 d\nP 12 2 dlnPn d In P 2 2 1 _ ^ ^ ^ 

f/i2 di?i2 -ffn — H22 dHu dH 22 

This gives three differential equations for three independent functions Pij(Hij). The 
solutions are Gaussian and have the product 

P(H) = Cexp[-A(H^ + H 2 22 + 2H 2 l2 ) - B(H n + H 22 )]- (2.14) 

The constant B can be eliminated by choosing an appropriate zero of energy. Then 
P(H) can be written as 

P(H) = Cexp(-ATr{H 2 }). (2.15) 

In a similar way, we can describe the construction of Gaussian Hermitian matri- 
ces where the probability P(H) is invariant under the unitary transformation. To 
summarize [34, 9], different ensembles of random matrices follow from demanding (1) 
invariance of P{H) under the different possible groups of canonical transformations 
and (2) complete statistical non-correlations between all matrix elements. We now 
discuss briefly the circular ensembles 



2.1.3 Circular Ensembles 



The Gaussian ensembles are useful in studying the properties of the energy levels 
of the system. To model the dynamical behaviour of chaotic systems, study of 
the properties of random unitaries is useful. As mentioned above, the ensemble 
of random unitaries are also known as the "circular ensembles". The ensemble of 
random matrices used to describe the unitary evolution unrestricted by the time 
reversal symmetry is called the Circular Unitary Ensemble (CUE). Similarly, the 
ensemble of random matrices used to describe the unitary evolution having a time 
reversal symmetry are known as the Circular Orthogonal Ensemble (COE). 
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To discuss the properties of the eigenvalues of random unitary matrices, it will be 
useful to look at an example of a 2 x 2 random unitary matrix from the orthogonal 
ensemble [10]. A general symmetric unitary matrix can be expressed as: 



Te ir > 



Re 1 "! 
Te ir > 



and T + R = 1. The phases 7 and n are assumed to be independent and uniformly 
distributed between and 2ir. The eigenvalues of the matrix lie on the unit circle in 
the complex plane and are given by: 



Ai 2 = ±e i{v±(l,) = e iai ' 



(2.16) 



where sin</> = V-Rsin(7 — T)). We see that the difference between the phases, a = 
a± — ct2 = vr — 20, depends only on the the difference 7 — 77. The probability 
distribution function of 7 — rj is a constant given by P(7 — if) 
distribution function of the phase difference is given by: 



1=-. Thus, the 



P(a) 



1 

2^ 



9(7 - rj) 



da 



sin(a/2) 



4tt v / j R-cos 2 («/2) 



(2.17) 



We see that as the phase differences become smaller, the probability distribution 
function is suppressed. This is commonly known as the repulsion of eigenvalues for 
random unitary matrices. Similar repulsion occurs for matrices in higher dimensions. 
Level repulsion also occurs for random unitary matrices invariant under the unitary 
transformation [10]. Properties of eigenvalues like the level repulsion, as predicted 
by random matrix theory, agree remarkably well with the Hamitonians and unitaries 
describing a quantized chaotic map. 
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2.2 Shannon and Fisher Information 

In information theory, entropy is a measure of the uncertainty associated with a 
random variable. In this context, the term usually refers to the Shannon entropy, 
which quantifies the expected value of the information contained in a message, usually 
in units such as bits. In this context, a "message" means a specific realization of 
the random variable. In classical information theory, the Shannon entropy, H is 
associated to a probability distribution in the following way: 



where {p(l),p(2), ...,p(i)...} are the probability values taken by some random vari- 
able. The Shannon entropy is a measure of the amount of information and is related 
to the physical resources required to solve certain information processing tasks. For 
example, Shannon's noiseless channel coding theorem [15] says the following - Sup- 
pose Xi is an i.i.d information source with entropy rate H(X) (bits per symbol). 
Suppose R > H(X). Then there exists a reliable compression scheme of rate R for 
the source. Conversely, if R < H(X) then any compression scheme will not be re- 
liable. Thus, Shannon's entropy represents an absolute limit on the best possible 
lossless compression of any communication. This operational interpretation of Shan- 
non entropy in terms of data compression is the cornerstone of classical information 
theory. 

The Fisher information is a way of measuring the amount of information that 
an observable random variable X carries about an unknown parameter upon which 
the probability of X depends. Thus, the Fisher information tells us how well we can 
estimate a parameter given a sequence of data. 

We start with a few definitions taken from [32]. Let {f(x; 9)}, 9 G G, denote an 
indexed family of probability densities, f(x; 9) > 0, J f(x; 6)dx = 1 for all 9. Here 9 
is the parameter set. 




(2.18) 
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Definition: An estimator for 9 for sample size n is a function T, whose input is 
the sequence of n data values and the output is an estimate, 9 est , of the parameter 
9. More formally, T : F n — > 0. The estimator gives us the approximate value of 
the parameter based upon the data X. It is desirable to have an idea of how good 
the estimator is in estimating the value of the parameter. Here, as we shall see, the 
Fisher information plays a crucial role here. 

The Fisher information is defined as 

J e = E e [-^\nf(x;9)] 2 . (2.19) 

The subscript 9 means that the expectation is with respect to the density /(.;#). 
The Fisher information gives a lower bound on the error in estimating 9 from the 
data. We state this result, which shows the significance of the Fisher information, 
without proof [32]. 

Theorem : The mean squared error of any unbiased estimator T(X) of the param- 
eter 9 is lower bounded by the reciprocal of the Fisher information, i.e., 

var(T) > ^. (2.20) 

This is known as the Cramer- Rao bound [35]. The bound states that the variance of 
any unbiased estimator is at least as high as the inverse of the Fisher information. 
An unbiased estimator which achieves this lower bound is said to be efficient. Hence, 
the Fisher information gives us the idea of the goodness of the estimator. 



2.3 Entanglement 

Entanglement is a uniquely quantum phenomena. It is the property of composite 
quantum mechanical systems that describes the non-separability of the constituents. 
Suppose, for example, Alice has an electron in spin up state, | f), while Bob has a 
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an electron in the spin down state | 1). Then we say that the joint state of Alice's 
and Bob's electron is, 

lt)®||> = |tl>- (2-21) 

These are called tensor product states. The Hilbert space that encompasses tensor 
product states is, however, not restricted by these states. For example the state, 
, known as the EPR pair (after Einstein, Podolsky and Rosen) is a per- 
fectly valid physical state. Entanglement tries to capture this inability to express a 
quantum state as a tensor product. Thus, entanglement can be regarded as a joint 
property of the system. Entanglement cannot be increased under local operations 
and classical communication and therefore has a non-local character to it. 

Quantum entanglement has shed light on some of the most fundamental aspects of 
quantum mechanics. For example, entangled states like the "EPR" pair can be used 
to show that "locality" (the assumption that measurements made at a location A 
cannot influence the measurement outcomes at another location B that is "space-like" 
separated from A) and "realism" (physical properties have definite values which exist 
independent of measurement) are incompatible with each other. The assumptions 
of "locality" and "realism" are together known as local realism. John Bell [36], in 
1964, demonstrated that if there is a realistic interpretation of quantum theory, in 
the sense defined above, it has to be nonlocal. Nature is not locally realistic! 

Not only does entanglement teach us important lessons in the foundations of 
quantum mechanics, it enables certain tasks which either cannot be performed clas- 
sically or would otherwise require far more resources. Entanglement has a crucial 
role in quantum teleportation [37], quantum dense coding [38], cryptography with 
the Bell theorem [40], and quantum computation using pure quantum states [41]. 

Now the question arises: how do we quantify the entanglement between two 
quantum systems? For this we look at an important result known as the Schmidt 
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decomposition in the next section. 

2.3.1 Quantifying Entanglement 

We consider entanglement of pure states of the bipartite system. Entanglement is 
then uniquely determined by the coefficients in the Schmidt decomposition of the 
joint state of the system, 



where A» are the eigenvalues of the reduced density matrix of either subsystem, 
and the Schmidt basis vectors {\ui)j, are their respective eigenvectors. The 

entanglement E is the Shannon entropy of the squares of the Schmidt coefficients, 



Unfortunately, quantification of entanglement for mixed states is not as simple. 
Extending the definition of separability from the pure case, we say that a mixed state 
is separable if it can be written as 



where p A and p B are themselves states on the subsystems A and B respectively. 
Therefore, a state is separable if it is a probability distribution over uncorrelated 
states, or product states. p A and p B can be taken to be pure states without any 
loss in generality. A state is then said to be entangled if it is not separable. As 
mentioned above, in general, finding out whether or not a mixed state is entangled 
is considered difficult. In general, the problem of determining whether a bipartite 
system in a mixed quantum state is entangled or not is NP hard [45]. For a detailed 
review of mixed state and multi-partite entanglement we refer the reader to [42]. 




(2.22) 




(2.23) 




(2.24) 
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The theory of entanglement has been very useful in helping us understand certain 
features of quantum mechanics. Entanglement also serves as a uniquely quantum re- 
source which can be harnessed to accomplish various quantum information processing 
tasks. This is the focus of our next section. 

2.4 The quantum advantage behind quantum in- 
formation processing 

The crowning achievement of quantum information science has been the discovery of 
algorithms and communication protocols based on the laws of quantum mechanics 
that accomplish tasks much more efficiently than the known classical methods. For 
example the field of quantum computation has supplied us with quantum algorithms 
for integer factorization [43, 116] and for searching an unsorted database [117] that 
offer significant improvements over the best known classical algorithms. 

Characterizing the resources behind the enhancements and speedups provided 
by quantum mechanics over best known classical procedures is one of the most fun- 
damental questions in quantum information science. Quantum entanglement [! '] is 
generally seen to be the key resource that gives quantum information processors their 
power. There are, however, quantum processes which provide an exponential advan- 
tage in the presence of little or no entanglement [98, 99]. In the realm of mixed-state 
quantum computation, for example, quantum discord [2, 50] has been proposed as a 
resource [3] and there has been progress in this direction in the last few years [100]. 
It has also been shown to be a resource in quantum state discrimination [101, 102] 
and quantum locking [103]. 

Quantum mechanics also provides an advantage in communication. For exam- 
ple, sending secret messages between two parties in the presence of adversaries in a 
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protocol known as quantum cryptography [40, 44] provides security well beyond any 
known classical protocol. The advantage of quantum cryptography is that it enables 
certain cryptographic tasks that are impossible using only classical communication. 
In particular, quantum mechanics guarantees that measuring quantum data disturbs 
that data; this can be used to detect an adversary's interference with a message. 

Quantum cryptography is not the only setting where we utilize the features of 
quantum mechanics to get an advantage over best known classical communication 
protocols. Research in quantum information science has led to a discovery of new 
protocols and tasks which are classically impossible. The most relevant example is 
that of quantum teleportation [37]. Quantum teleportation, or entanglement-assisted 
teleportation, is the process by which an unknown quantum state can be transmit- 
ted from one location to another, without the state being transmitted through the 
intervening space. 

We first introduce the notation typically used in quantum communication pro- 
tocols: [q q] represents one qubit of communication between two parties and [qq] 
represents one shared ebit (one maximally entangled two qubit state) between two 
parties. Similarly, [c — > c] represents one classical bit of communication between the 
parties. Expressing teleportation as a resource inequality [ )] we get, 



This inequality says that one can employ a shared ebit and 2 bits of classical com- 
munication to communicate a single unknown quantum bit. To communicate an 
unknown quantum bit by a classical procedure will take exponentially large amount 
of resources as compared to teleportation. 

Another protocol which shows the advantage of quantum communications is the 
superdense coding [38]. Expressing it as a resource inequality [30], 



[qq} + 2[c^c] y [q q). 




[qq] + [q->q] h 2[c ^ c], 



(2.26) 
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showing that one can employ a shared ebit and a single bit of quantum commu- 
nication to communicate 2 bits of classical information. Thus, complete classical 
information about two particles can be sent by direct manipulation of just one parti- 
cle by the sender. This protocol starts with Alice and Bob sharing a Bell state. For 
example, they share the state Jll^ti). Alice then performs one of the following op- 
erations on her half of the Bell state - either the identity operation, J, or applies one 
of the three operators, i.e. a x ,i<j y , or a z . The type of operation that Alice performs 
thus represents two bits of information which she wants to communicate with Bob. 
This she does by sending over her half of the Bell state to Bob after performing the 
operation. Depending on Alice's operation (one of J, a x ,i<r y , or a z ), Bob gets one 
of the four states: ltt) ^ u) , Jitt+M ln) ^ lt \ or ln) ^ u) . These states form what is 
known as the Bell basis, which is an orthonormal basis. Since orthonormal quantum 
states can be distinguished by a suitable quantum measurement, Bob comes to know 
of the joint state he shared with Alice and hence the operation performed by Alice 
at her end. Thus, Alice succeeds in communicating two bits of classical information 
to Bob using a shared ebit and a single bit of quantum communication. 

In both the protocols mentioned above, it is worth noting that it was the pre- 
existing entanglement that made them possible. Without entanglement, it is impos- 
sible to achieve both the teleportation and the super-dense coding. 

Thus we have seen a spectrum of quantum information processing protocols, 
from quantum computation to communications tasks and cryptography. Certain 
fundamental properties of quantum mechanics are exploited in accomplishing these 
tasks [ ]. These include: 

1. Superposition: A quantum state can exist in an arbitrary complex linear com- 
bination of classical logic states. A classic example is the "cat state" , named 
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after Schrodinger. 

\Alive) + \Dead) 

Wcat) = ^ (2.27) 

Both of the logic states evolve in "parallel" according to a given unitary evo- 
lution. 

2. Interference: The quantum wave function undergoes interference, just like the 
classical waves. The quantum state can be in a superposition of different 
alternatives/paths, each of which is characterized by a complex probability 
amplitude, that has both real and imaginary components in general. If these 
paths are in principle indistinguishable, i.e. there is no information whatsoever 
in the universe that can help us distinguish between them, then in that case 
the corresponding probability amplitudes add up causing interference. This 
feature is utilized in quantum computers, where different paths are explored 
in parallel in search of the solution and the probability amplitude of the path 
leading to the right solution gradually builds up. 

3. Entanglement: As mentioned earlier, protocols like teleportation, superdense 
coding and pure state quantum computation are made possible due to entan- 
glement. 

4. Non-distinguishability of quantum states: Non-orthogonal quantum states can- 
not be unambiguously distinguished. Moreover, obtaining information about 
an unknown quantum state can cause disturbance and actually change it. This 
feature is exploited in designing cryptographic protocols. 

As we have seen, the theory of entanglement has been very helpful in understand- 
ing quantum mechanics from a fundamental point of view. In addition, entanglement 
has been shown to be a crucial resource in quantum information processing. How- 
ever, this is not the whole story. In recent years, there has been some progress in 
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quantifying the quantum character of composite quantum systems that goes beyond 
entanglement. This is the focus of our next section. 

2.5 Quantum Discord: A quantum resource? 

Quantum discord aims at generalizing the notion of quantum correlations in a quan- 
tum state, beyond entanglement [2, 50, 51]. It aims to capture all the nonclassical 
correlations in a quantum system. Quantum measurements disturb a quantum sys- 
tem in a way that is unique to quantum theory. Quantum correlations in a bipartite 
system are precisely the ones that are destroyed by such disturbances. As we discuss 
below, this feature of quantum systems can be used to quantify the amount of purely 
quantum correlations present in a bipartite quantum system. 

Quantum mutual information is generally taken to be the measure of total cor- 
relations, classical and quantum, in a quantum state. For two systems, A and B, 
it is defined as I(A : B) = H(A) + H(B) - H(A, B). Here #(-) denotes the von 
Neumann entropy of the appropriate distribution [ 9]. For a classical probability 
distribution, Bayes' rule leads to an equivalent definition of the mutual information 
as I(A : B) = H(A) — H(A\B). This motivates a definition of classical correlation in 
a quantum state. 

Suppose Alice and Bob share a quantum state pab £ Ha ®Hb- If Bob performs 
a measurement specified by the POVM set {IT}, the resulting state is given by the 
shared ensemble {pi,PA\i}, where 

p A \i = Tr B (UipAB)/Pi, Pi = Tr a ,b (Jli Pab)- 

A quantum analogue of the conditional entropy can then be defined as S{jii}(A\B) = 
^2iPi.S(pA\i), and an alternative version of the quantum mutual information can 
now be defined as J{u 1 }{.Pab) = S(pa) — S{iii}(A\B), where S(-) denotes the von 



30 



Chapter 2. Background 



Neumann entropy of the relevant state. The above quantity depends on the chosen 
set of measurements {IT}. To capture all the classical correlations present in Pab, we 
maximize ^{^(pab) over all {IT}, arriving at a measurement independent quantity 



J( PAB ) = nux{S(p A ) ~ S { n l} (A\B)). (2.28) 



Then, quantum discord is defined as [ 



V(p AB ) = I (p AB ) - J (Pab) (2.29) 
= S(p B ) - S(pab) +min S {Ui} (A\B) 



2.5.1 Operational interpretation of quantum discord 

The work in this dissertation seeks to answer very fundamental question about quan- 
tum discord: Is discord just a mathematical construct or does it have a definite 
physical role in information processing? We provide a link between quantum discord 
and an actual physical task involving communication between two parties [4] . 

The key insight to our findings is that quantum measurements and environmental 
decoherence disturb a quantum system in a way that is unique to quantum theory. 
Quantum correlations in a bipartite system are precisely the ones that are destroyed 
by such disturbances, and therefore certain quantum communication protocols be- 
come overloaded by an amount exactly equal to quantum discord. More specifically, 
we will show in Chapter 5 that discord is the markup in the cost of quantum com- 
munication in the process of quantum state merging [114], if the system undergoes 
measurement and/or decoherence. 

At this point, a natural question arises: Is there a role of quantum discord in 
quantum information theory as a whole beyond the lossy state merging protocol? 



31 



Chapter 2. Background 

We answer this question with a yes. We first observe that quantum state merging 
protocol is a derivative of the more general Fully Quantum Slepian Wolf (FQSW) 
protocol [114, 3 ] and the closely related "mother" protocol. The mother protocol 
is essentially a unification of all unidirectional, bipartite and memoryless quantum 
communication protocols like quantum teleportation, superdense coding and entan- 
glement distillation. A link between discord and state merging does indeed suggest 
a link between discord and the mother protocol and hence a possible role of discord 
in all bipartite, unidirectional and memoryless quantum communication protocols. 

Before exploring the role of discord in quantum information theory, it is therefore 
necessary to review the mother protocol and other related concepts. 

2.6 The mother protocol and quantum informa- 
tion's family tree 

So far, we have been discussing quantum communication tasks like teleportation, 
super-dense coding, quantum state merging as independent entities. Is there a uni- 
fied way of looking at various quantum communication tasks? Do various quantum 
communication protocols described above have a common origin. Interestingly, the 
answer is yes. It was shown in [33] that essentially all unidirectional, bipartite and 
memoryless quantum communication protocols are actually siblings originating from 
one "mother" . The mother protocol can be seen to provide a hierarchical structure to 
the family of quantum protocols. We will describe the improved version of the mother 
protocol, which is the fully quantum Slepian- Wolf (FQSW) protocol. Throughout 
the thesis, we will use the names "mother" and FQSW to describe essentially this 
same protocol. 
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2.6.1 The mother protocol 

At the start of the mother protocol, we have n copies of a quantum state that can 
be expressed as the tensor power \(v() R )® n ) . Alice holds the A shares and Bob the 
B shares. The reference system R is "purification" of the AB system (which might 
be described by a mixed state) and does not actively participate in the protocol. 
The mother protocol can be viewed as an entanglement distillation between A and 
B when the only type of communication permitted is the ability to send qubits from 
Alice to Bob. The transformation can be expressed in the resource inequality as 



0j AB ) + l -I{A : R)[q -+q}> l -I{A : B)[qq}. (2.30) 

Here, \^ AB ) refers to the state shared between Alice and Bob whose purification 
is the state \t/j ABR ). The above inequality states that n copies of the state |?/> AB ) 
can be converted to \l{A : B) EPR pairs per copy, provided Alice is allowed to 
communicate with Bob by sending him qubits at the rate \I{A : R) per copy. 

The mother protocol has a stronger version known as the FQSW protocol. The 
FQSW protocol not only enables the two parties, Alice (A) and Bob (B), to distill 
\l{A : B) EPR pairs per copy, in addition Alice can "merge" her state with Bob. 
By the state merging task, we mean that Alice is able to successfully transfer her 
entanglement with the reference system R to Bob as shown in Fig 2.1. 

Writing the FQSW in terms of a resource inequality 

I \ 
{^ AB } + -I(A :R)[q-+q]> -I (A : B) [gg]+State Merging between A and B (2.31) 

In a more rigorous mathematical notation, we write the above as 

(U S ^ AB : ^ S ) + \l{A :R)[q^q}> \l(A : B)[qq] + (id s ^ B : ^ s ). (2.32) 
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Figure 2.1: a) The starting point of the FQSW protocol. A pure tripartite entangled 
state \(ip ABR ^® n } shared between Alice and Bob. R is the reference system purifying 
AB system and does not participate actively in the protocol, b) After execution of 
the protocol. Alice's portion of the state has been transferred to Bob, and now Bob 
and R hold the state \(if> RB ) 9n ) where B is a register held with B. 



It is important to explain the terminology used in the above inequality. When we 
have a noisy resource like a mixed state, ip s , or a noisy channel, it is inserted between 
a "()". Thus a mixed state is represented by (V^), and a noisy channel by (A/"). A 
channel is a relative resource meaning that the protocol only works 

provided the input to the channel is the state tf) s . On the LHS, hi takes the state 
ip s and distributes it to Alice and Bob. On the RHS, the symbol "zd" is an identity 
channel taking the state ip s to Bob alone. The state ?/> 5 on the left-hand side of 
the inequality is distributed to Alice and Bob, while on the right-hand side, that 
same state is given to Bob alone. This inequality states that starting from the state 
\(ip ABR )® n ) , and using \I(A : R) bits of quantum communication from Alice to Bob, 
they can distill hi (A : B) EPR pairs per copy, and in addition Alice can accomplish 
merging her state with Bob. Figure 2.1 [33] shows the action of the protocol. 

In the state merging task, as described above, Alice is able to successfully trans- 
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fer her entanglement with the reference system R to Bob. This means that Alice 
transfers her portion of the state to Bob. In other words, they manage to create the 
state 10^)®"), where B is a register held with B and \(ip R§ ) m ) = \{^ ABR )® n ) in 
the limit n — > oo. 

Thus the FQSW protocol accomplishes state merging as well as entanglement 
distillation. The FQSW protocol is valid asymptotically in the limit of a large number 
of copies and this is denoted by the symbol >. 
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Dynamically Generated 
Entanglement and Chaos 



3.1 Introduction 



The connections between complexity, nonlinear dynamics, ergodicity, and entropy 
production, have long been at the heart of the foundations of statistical physics. A 
central goal of "quantum chaos" has been to extend this foundation to the quan- 
tum world. Classic works on the subject including level statistics [46], properties of 
Wigner functions [4 ], and quantum scars in ergodic phase spaces [ 5] have tended 
to focus on the properties of wave mechanics, e.g. the dynamics of single particle 
billiards [ ] (also seen in the properties of classical waves, e.g. microwave cavities 
[49]). More recently, the tensor product structure of quantum mechanics, essential 
for understanding systems with multiple degrees of freedom, has come to the fore. 
In that context, one is naturally led to consider how the dynamical generation of 
entanglement between quantum subsystems is connected with the chaotic dynamics 
of coupled classical degrees of freedom. Such studies address fundamental issues of 
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complexity in quantum systems and are potentially applicable in quantum informa- 
tion processing, where entanglement is considered to be an essential resource. The 
connection between chaos in the classical description of Hamiltonian dynamics and 
entanglement in the quantum description has been the subject of extensive study 
over the last decade. The original motivation of Zurek and Paz was to address the 
quantum-to-classical transition [ ]. By conjecturing that chaotic systems decohere 
exponentially fast through their entanglement with the environment, they hoped to 
resolve a paradox in which a macroscopic system would exhibit the effects of quantum 
coherence on a time scale logarithmic in h. 

Work quickly following this turned to studies of the coupling of just two degrees 
of freedom, rather than system-environment coupling, as entanglement is most easily 
quantified for bipartite systems [56]. In most cases, workers have considered systems 
described by a total Hamiltonian of the form 

iWt) = Ht{t) + H 2 (t) + H int {t), (3.1) 

where Hi(t) and H2(t) can exhibit chaos in the classical description of the dynamics 
and H int (t) couples the two degrees of freedom. The most well-studied example has 
been two coupled kicked tops (a standard paradigm of quantum chaos [9]). Two sep- 
arate questions have been addressed: (i) How does the rate of dynamical generation 
of entanglement correlate with the chaos in the subsystems 1 and 2? (ii) How does 
the entanglement content of the state, either in the eigenstates, or in the state that 
is dynamically generated in quasi-steady state, correlate with this chaos? Miller and 
Sarkar [ 7] were the first to study question (i) for this system, and through numer- 
ical studies, correlated the rate of generation of entanglement with the Lyapunov 
exponents associated with the mean positions of quantum wavepackets localized in 
a mixed phase space (weak chaos). This behavior was not found to be universal 
[61, 59, 60, 58], but depended strongly on the degree of the chaos within subsystems 
when compared to the size of coupling between them. In a seeming paradox, for 



37 



Chapter 3. Dynamically Generated Entanglement and Chaos 

strong chaos within the tops, the rate of entanglement generation decreased with 
coupling strength, asymptoting to a constant value. Moreover, in a systematic study 
of the entangling power of the coupled kicked top system [62], Demkowicz-Dobrzaski 
and Kus found additional anomalies, including a regime in which each top was de- 
scribed by highly regular dynamics, but exhibited the highest rate of generation of 
entanglement when compared with conditions where each top is highly chaotic. 

In relation to question (ii), Bandyopadhyay and Lakshminarayan [63, 58] ex- 
plored the amount of entanglement that is associated with coupled kicked tops, with 
particular emphasis on the entanglement of the Floquet eigenstates [64]. The en- 
tanglement of these eigenstates saturated to a value below the maximum possible 
value in a way that depended only on the Hilbert space dimension, not the chaoticity 
parameter. The same was true of the dynamically generated entanglement. (The 
relationship between the entanglement in the eigenstates and the dynamically gener- 
ated entanglement is subtle [62]; we'll return to this point later). This work gave the 
first indication that the entanglement generated by the coupled tops was statistical 
in nature, and related to the theory of random states in Hilbert space. Using ran- 
dom matrix theory [34, 9] they were able to determine the statistics of the Schmidt 
coefficients of a random bipartite pure state, and thus were able to predict the satu- 
ration value of the entanglement for the Floquet eigenstates. An extended analysis 
of the statistical properties of the Schmidt vectors of random states was carried out 
by Znidaric [65]. In other related work, dynamical generation of entanglement by 
chaotic maps and its relation to random matrices was also explored by Gorin and 
Seligman [66] as a way of modeling decoherence, by Scott and Caves [21] and Abreu 
and Vallejos [67] as a way of comparing different quantizations of the Baker's map, 
and by Viola and coworkers [68] as a means of quantifying complexity in quantum 
systems and its relationship to generalized entanglement. 

In another approach to question (ii), Ghose and Sanders have shown that there 
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are signatures of chaos in the entanglement dynamically generated by a single kicked 
top when the large angular momentum is thought of as a collection of symmetri- 
cally coupled qubits [69, 70]. They showed strong correlation between the classical 
Poincare surface of section for a mixed phase space, and a contour plot of the dy- 
namically generated entanglement as a function of the initial position of a localized 
coherent state. Using the Floquet spectrum, they also explained the initial rise time 
and power spectrum in the entanglement history. 

While many of the elements connecting chaos and entanglement have been ex- 
plored with a variety of successful numerical and analytic predictions, in some cases 
the key relations have been obscured. In particular, in studies of systems of the form 
in Eq. 3.1, entanglement is correlated with the chaos in the individual subsystems. 
But, entanglement arises from coupling between subsystems, and is a global property 
of the state. Likewise, chaos can also arise through the coupling of degrees of freedom 
when the overall dynamics are not integrable. For this reason, we believe the key 
relations are best understood by correlating entanglement with chaos in the joint 
system (i.e., chaos in H tota \), rather than chaos in the subsystems that one would 
see in the absence of coupling. To do so, it is most natural to consider systems in 
which chaos and entanglement arise from the same mechanism - the physical cou- 
pling between subsystems. Moreover, by considering a total system that is chaotic 
only when the two parts are coupled, we focus on a classical phase space that de- 
scribes the global system rather than a subsystem, and there is no ambiguity about 
the nature of the joint dynamics. For this case, the distinction between weak and 
strong coupling cannot be made independently of weak and strong chaos, thereby 
sharpening our focus on the key relationships. 

To address these issues, we consider a model system of kicked coupled-tops, rather 
then coupled kicked-tops, described in detail in this chapter. This system is moti- 
vated by its connection to possible experimental realizations, our ability to easily 
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visualize the classical phase space, and to analyze the Floquet map. We use this 
system as a forum to explore question (ii) above - how is chaos in the classical 
dynamics of the joint system correlated with the long-time averaged dynamically 
generated entanglement? 

The basic thesis of the work described in this chapter is as follows. Chaos can arise 
in classical dynamics when there are insufficient symmetries (integrals of motion) for a 
given number of degrees of freedom. In the quantum analog, insufficient symmetries 
lead to the random matrix conjecture - systems with global classical chaos have 
eigenvectors and eigenvalues that are statistically predicted by ensembles of random 
matrices [11, 9]. Moreover, where classical chaos leads to ergodic dynamics and 
the generation of "random" coarse-grained distributions on phase space, for times 
short compared to the Heisenberg time, but long compared to transient behavior, the 
quantum chaotic map generates a state with many properties that are statistically 
predicted by a random state in Hilbert space, picked according to the appropriate 
Haar measure [21]. The dynamically generated entanglement is then that of a random 
state (by this measure) in the relevant Hilbert state. These predictions can be 
extended to mixed phase spaces with regular islands immersed in a chaotic sea. 
With the help of Percival's conjecture [71] that divides eigenstates into chaotic and 
regular classes, we can find the entanglement of a random state in a chaotic subspace 
and thus predict entanglement generation in a mixed phase space for chaotic initial 
conditions. Whereas in the globally chaotic case we can derive analytic results, for 
the mixed phase spaces we are relegated to numerical predictions, which nonetheless 
verify the connection between entanglement generation in chaotic dynamics and the 
creation of pseudo-random states in Hilbert space. 

The remainder of this chapter is organized as follows. In Sec. 3.2 we introduce our 
model of kicked coupled-tops, studying the classical and quantum features. Section 
3.3, the heart of this work, studies the entanglement in our system. We perform 
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numerical calculations of the entanglement of the system's eigenstates, the long- 
time averaged entanglement generated by the Floquet map, and its relationship to 
the classical phase space. We then explain these results in terms of the properties of 
random states in Hilbert space. Reviewing the essential ideas, we derive new analytic 
expressions for the typical entanglement of a random state when we are restricted 
to a subspace of the full tensor product space. This is of relevance here given the 
symmetries of the system. We also pay particular attention to the subtle distinctions 
between the eigenstates of random matrices and the random states generated from 
initially localized wavepackets. In doing so we clarify previous works and make 
accurate predictions, especially for global chaos, but also extended to a more general 
mixed phase space scenario. Our results are discussed and summarized in Sec. 3.4. 



3.2.1 Quantum and classical descriptions 

We consider a bipartite system composed of two spins, I and J, isotropically coupled 
in a Heisenberg interaction, and subject to periodic kicks that act only on spin 
J. Choosing the direction of the kicks to be about the z-axis, the system evolves 
according to the Hamiltonian, 



Here A gives the strength of the isotropic coupling, B the strength of the kicking, 
and r is the kicking period. Such a Hamiltonian describes the hyperfme interaction 
between nuclear spin I and total electron angular momentum J, with a magnetic field 
that has negligible effect on the nucleus. While this realization cannot reach deep 
into the semiclassical regime, for large atoms, with heavy nuclei and a large number 



3.2 Kicked Coupled- Tops 



oo 




(3.2) 



n=— oo 
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of electrons in the valance shell, one can explore nontrival mesoscopic regimes. The 
true semiclassical limit can potentially be attained in an atom-photon system where 
I is the collective spin of an atomic ensemble coupled to the Stokes vector J of 
a quantized electromagnetic field [' ]. We will not consider here the feasibility of 
experimental realizations, instead focusing on the foundational theory. 

Choosing the external field to act in delta kicks allows us to express the Floquet 
map (transformation after one period) in a simple form of sequential rotations, 

U T = e -*M e -W. = e -iaFV2 e -i0J^ (33) 

where a and /3 are related to A and B in terms of the kicking period. The Floquet 
map describes the evolution of the system for one period. In the second form, we 
have expressed the rotation in terms of the total angular momentum F = I + J 
and neglected irrelevant overall phases. We can thus interpret the dynamics as 
alternating a rotation of J about a space fixed z-axis by angle /3, followed by a 
procession of I and J about F by an angle a|F|, as shown in Fig (3.1). Such a 
simple transformation nonetheless leads to complex dynamics, including chaos in 
the classical limit as discussed below. From the quantum perspective, since the 
two rotations don't commute, there are insufficient symmetries to specify Floquet 
eigenstates by a complete set of commuting operators; the system is not integrable. 
Note, however, that the system is invariant under an overall rotation around the 
z-axis, so F z is a conserved quantity (F 2 is not conserved). 

We treat the classical limit of quantum mechanical spin in the familiar way [9]. 
Each of our spins has three components, but a fixed magnitude, and thus their orien- 
tations can be specified by two variables. The z-component of a spin and the angle 
4>, denoting its orientation in the x-y plane, are canonically conjugate, and thus each 
spin constitutes one canonical degree of freedom. The classical dynamical map has 
the same physical action as described above in the quantum context - rotation of J 
by angle /3 followed by precession of I and J about F by angle a|F|. Here, the rota- 
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B A 




Figure 3.1: The dynamics of the kicked coupled tops can be viewed as an alternating 
sequence of rotations. The two spins I and J precess around the total angular 
momentum F by an angle a|F|, and the spin J is kicked around the space-fixed 
z-axis by (3. 



tions are implemented by 3 x 3 SO (3) matrices. The two spins, plus time-dependent 
Hamiltonian imply a five dimensional phase space. Since F z is conserved, the dy- 
namics is restricted to a four- dimensional hypersurface. As there are no additional 
constraints, the dynamics are not integrable and can exhibit chaos. Note, Eq. (3.2) 
is of the form of Eq. 3.1, with H2 = 0, but where chaos is only seen in the coupled 
dynamics, not the dynamics of the of Hi alone. 



To visualize the dynamics, we rewrite our system in terms of a new set of variables, 
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(F z , = 0/ + 0j) and (5F Z = I z - J z , 50 = 0/ - 0j), 

F z - 5F Z F Z + 5F Z 
Jz = 2 ' Iz = 2 ' ^ ^ 

I • J = I z J z + \J I 2 — I 2 \J ' J 2 — J 2 (sin 0/ sin 0j + cos 0/ cos 0j) 



(3.4b) 

Because F z is a conserved quantity, does not appear in our Hamiltonian. It is a 
cyclic coordinate, and thus we can ignore it without losing any information about the 
further evolution of the remaining variables. Neither do we require to determine 
the Lyapunov exponent of a chaotic system. Thus, we need only consider the two 
difference variables, (5F Z , 50), and time, taking us from a four to a three dimensional 
hypersurface. This allows us to visualize our system using a Poincare surface of 
section as a stroboscopic plot. We restrict our attention here to F z = as this also 
leads to the largest subspace in the associated quantum problem. The reduction 
of our system to essentially one degree of freedom is not generic, but simplifies the 
analysis without sacrificing our ability to study the essential relations between chaos 
and entanglement. 



cos 
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Figure 3.2: Poincare surface of section for the coupled kicked tops, with F z 
Regular phase motion: a — 1/2, /3 — tt/2, (b) Mixed phase space: a = 3/2, (5 
(c) Global chaos: a = 6,/3 = 7r/2. 
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The classical equations of motion depend on the ratio |I|/|J|. We focus here on 
equal spin magnitudes and fix F z = 0. Thus, without loss of generality, since the 
SO (3) rotation matrices of classical dynamics are independent of spin magnitude, we 
take the spin vectors to be unit vectors. The basic structure of the phase space can be 
understood as follows. When the coupling is removed, our system has fixed points 
at the northern and southern "poles". As the chaoticity parameter is turned up, 
chaos first forms around the unstable "north pole" while regular behavior persists 
around the stable "south pole". Further fixed points appear in the usual manner 
as bifurcations occur with increase of the chaoticity parameter. Figure (3.2) shows 
three different regimes of classical dynamics. With the parameters a = 1/2, (3 = 7r/2 
(Fig. (3.2a), the dynamics are highly regular, with negligible stochastic motion. 
When a = 3/2, /3 = n/2 (Fig. (3.2b), we see a mixed space with chaotic and regular 
regions of comparable size. The parameters a = 6, /3 = tt /2 (Fig. (3.2c), give a 
completely chaotic phase space. 

We want to choose our quantum Hamiltonian so that we will recover our classical 
dynamics in the large spin limit. We would like to be able to vary the size of our 
spins, but we will keep the pair equal to each other in magnitude, I = J. Since the 
SU(2) rotation matrices depend on the spin magnitude, we must scale the Floquet 
operator. By substituting a — >a = a/Jwe obtain the same Heisenberg equations 
of motion as the classical equations for equal magnitude spins. 

3.2.2 Quantum chaology 

In order to understand the dynamical generation of entanglement, we need to estab- 
lish some basic understanding of the eigenstates of the system and their relationship 
to the classical dynamics. As our system is time periodic, the states of interest are 
the eigenstates of the Floquet operator, Eq. (3.3). It is useful to consider both the 
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coupled and uncoupled representations of angular momentum connected by the usual 
Clebsch-Gordan expansion, 

\F,M F )= ^ (I,mr,J,mj\F,M F )\I,mi)\J,mj) (3.5) 

mi,mj 

Conservation of F z implies that the operator is block diagonal for all states defined 
by quantum number Mp. The largest block, Mp = 0, has dimension 2 J + 1 as F 
varies from to 2 J. Using the uncoupled representation, denoting the product state 
by the single quantum number mj = —mj, the matrix, 

(m'j\U T \mj) = ^^("^iir^+^-O (3.6) 

F 

(F, 0\I, -m'j- J, m'j) (F, 0\I, -my, J, mj) 

can then be diagonalized to yield the Floquet eigenstates and eigenphases, 

\k) = 2cW|J,-mj)|J,mj); U T \k) = e^\k). (3.7) 

mj 

A central result of quantum chaos is the connection with the theory of random 
matrices [9]. In the limit of large Hilbert space dimensions (small h), for parameters 
such that the classical description of the dynamics shows global chaos, the eigen- 
states and eigenvalues of the quantum dynamics have the statistical properties of an 
ensemble of random matrices. The appropriate ensemble depends on the properties 
of the quantum system under time-reversal [ ]. We thus seek to determine whether 
there exists an anti-unitary (time reversal) operator T that has the following action 
on the Floquet operator, 

TU T T~ l = U ] T = e i/3J *e i&I - J . (3.8) 

Analogous to the case of the single kicked top, we consider the generalized time 
reversal operation, 

T = e il3 ' h K, (3.9) 
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where K is complex conjugation in the uncoupled product representation. Since both 
I y and J y change sign under conjugation, while the x and z components do not, 

KJ Z K = J Z - KIJK = IJ. (3.10) 

It then follows that 

TU T T- 1 = (e^ Jz K) (e- i&LJ e- i(3Jz ) (Ke~ m *) (3.11) 
= e i/3Jz (e i&I ' 3 e^ Jz ) e~ ifiJz 
= e ll3Jz e i&I - J = Ul 

so the dynamics are time-reversal invariant. Moreover, T 2 = 1, so there is no 
Kramer's degeneracy. Given these facts, for parameters in which the classical dy- 
namics are globally chaotic, we expect the Floquet operator to have the statistical 
properties of a random matrix chosen from the circular orthogonal ensemble (COE). 

To further correlate the Floquet eigenstates with the classical phase space in the 
case of regular and mixed dynamics, it is useful to employ a Husimi representation. 
A spin coherent state has a minimum quantum uncertainty and is specified by polar 
orientation angles 9 and <fi on the sphere. In terms of the standard basis, a spin 
coherent state for a single spin is [73] 



\n) = Y — r J-n — rfr! d J >™>> ( 3 - 12 ) 

where \x = tan(#/2)e J9i . For our system, because the subspaces in which the eigen- 
states live are not described by an irreducible representation of angular momentum, 
there are no such minimum uncertainty states for the difference angles. Nonetheless, 
we obtain a useful set of states by projecting the product of spin coherent states 
associated with the two subsystems onto the subspace with a fixed value of F z (here 
F z — 0). The result of the projection is 

PoMM = E (-Yn ( wL d m) f l-m)j. (3-13) 

^\^JJ \J — m)\{J + m)\ 
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Classically, in projecting onto the surface of section with F z = 0, we take 8i + 8j = ir. 
Fixing this value in the quantum state one finds 

"l + sin 



i(<f>i-<Ps) 



(3.14) 



l-sin(^) 

The projected coherent state thus depends only on the difference of the angle vari- 
ables, and allows us to consider localized quantum states correlated with the classical 
phase space of interest. After normalizing, we arrive at an over-complete basis of 
states for the F z = subspace, parameterized by 58 and S(f> . The Husimi distribution 
of a state \ip) in this space, 



(3.15) 



then provides a visualization in phase space. 
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Figure 3.3: Husimi distributions of Floquet eigenstates associated with the parame- 
ters of a mixed phase space (Fig. 3.2b). (a,b,c) Regular eigenstates around different 
fixed points. (d,e,f ) Chaotic eigenstates, delocalized in the chaotic sea. For spin size, 
we have I = J = 150. 
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In order to explore the semiclassical limit, we choose I = J = 150, corresponding 
to a d = 301 dimensional Hilbert space in the F z = subspace, or an "effective H" 
of h e g = 1/301. A 301 dimensional Hilbert space is sufficiently large so that the 
minimum uncertainty spin coherent states can be effective considered as "points" 
on the phase space. As the size of the Hilbert space increases, we expect a better 
quantum-classical correspondence. Figure 3.3 shows the Husimi plots of a few of the 
eigenstates for a/ J = 3/2, (3 = ir/2, for which the classical phase space is mixed (Fig. 
3.2b). These plots exhibit the features expected according to Percival's conjecture. 
The states roughly divide into regular and irregular sets, with regular eigenstates 
concentrated on invariant tori around stable fixed points, resembling harmonic os- 
cillator eigenstates, and irregular "chaotic" states randomly distributed within the 
chaotic sea. 
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Figure 3.4: Scatter plot of (J z ) vs. Husimi Entropy, Sq, for the Floquet eigenstates 
associated with the mixed phase space (Fig. (3.2b)). Boxed regions (a), (b), and (c) 
correspond to regular states centered around fixed points. States in region (d) are 
considered "chaotic eigenstates". 
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Though Percival's conjecture is largely born out in numerical analyses, it is not 
strictly true (especially in the finite h limit), nor is there a strict procedure for 
filtering the regular from chaotic eigenstates except for very special systems [74]. We 
can, nonetheless, create an approximate filter. A useful measure for distinguishing 
states is the Shannon entropy of the Husimi distribution [75], 



where dfi is the measure on the phase space of difference angles on the sphere and Q 
is defined as in Eq. 3.15. We expect the states delocalized in the chaotic sea to have 
large entropy by this measure, while those states well-localized around fixed points 
have low entropy. This leaves some ambiguous situations, since highly excited states 
on regular tori also have high "Husimi entropy". To improve the filter, we follow a 
procedure suggested by Korsch and coworkers [76], which correlated the properties 
of the eigenstates to the classical phase space in order to distinguish the regular and 
irregular states for a nonlinear rotor. In Fig. 3.4 we plot the values of Sq and (J z )- 
The latter quantity correlates to the mean value of 56 in the semiclassical limit. We 
see four distinct features in this plot. Two lines of states with near constant ( J z ) but 
increasing Sq, boxed in Figs. (3.4a,b), correspond to the series of states localized 
around fixed points with increasing excitation (Figs 3.3a,b). The line of states with 
near constant Sq and increasing values of (J z ), boxed in Fig. 3.4c, correspond to 
the series of states localized around the stable "south pole" (Fig. 3.3c). Finally, 
the cluster of states with high values of both Sq and (J z ), boxed in Fig. 3.4d, 
correspond to the states delocalized in the chaotic sea that are concentrated near the 
original unstable fixed point at the "north pole" of the regular dynamics. There is 
no clean division between this cluster and states clearly localized on invariant tori. 
A qualitative examination, denoted in Fig 3.4, nonetheless gives us an indication 
of the chaotic subspace for these mixed dynamics. As the size of the Hilbert space 
increases, we expect a sharper distinction between regular and chaotic eigenstates. 
Such an identification is useful for giving quantitative prediction of the dynamically 




(3.16) 
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generated entanglement, as we discuss in the next section. 

3.3 Entanglement 

3.3.1 Calculating Entanglement 

We consider only pure states of the bipartite system. Entanglement is then uniquely 
determined by the coefficients in the Schmidt decomposition of the joint state of the 
system, 



where Aj are the eigenvalues of the reduced density matrix of either subsystem, 
and the Schmidt basis vectors are their respective eigenvectors. The 

entanglement E is the Shannon entropy of the Schmidt coefficients, 



Determination of the Schmidt decomposition is typically a nontrivial task, re- 
quiring partial trace and diagonalization of the reduced density operator [29]. The 
Schmidt basis will generally depend on the state \^)u- For the system at hand, we 
have a unique situation - within a subspace with a fixed value of E z , the uncoupled 
basis of angular momentum is the Schmidt basis, independent of the state, as seen, 
e.g., Eq. (3.7). Thus, for states within such subspaces, the entanglement is easily 
calculated as the Shannon entropy of the probability distribution of the state when 
expanded in the standard product basis. This not only simplifies calculations, but 
connects entanglement with the entropy of random states with respect to a fixed 
basis [77]. Thus, the choice of this particular metric of entanglement enables us to 




(3.17) 




(3.18) 
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compare the entanglement properties of a state to the properties of random states 
in Hilbert space. 

Throughout this section, we consider the F z = subspace, and take I = J = 150, 
corresponding to a Hilbert space of dimension d = 301. The maximum possible 
entanglement in this case is E max = logd m 5.71. 

3.3.2 Numerical Solutions 

The entanglement of the Floquet eigenstates is easily calculated based on the discus- 
sion above. Since the eigenstates reside in a subspace with fixed F z , the uncoupled 
representation of angular momentum is the Schmidt basis, and the entanglement 
between spins in a given eigenstate \k) is the Shannon entropy of the probability 
distribution of the expansion Amj = \cm\\ 2 from Eq. (3.7). Figure 3.5 shows a list 
plot of this entanglement for a mixed phase space (as shown in Fig. (3.2b)) and a 
completely chaotic space (as shown in Fig. (3.2c)). In the latter case, the entan- 
glement values are clustered around the value expected from random matrix theory, 
discussed below. 
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Figure 3.5: Entanglement of the Floquet eigenstates. (a) Map corresponding to a 
mixed phase space: a = 3/2,(3 = ir/2. (b) Map corresponding to global chaos: 
a = 6, (3 = 7r/2. The solid line gives the value expected from random matrix theory, 
Eq. (3.25). 
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Our main interest is to study the dynamically generated entanglement and its 
correlation with the classical phase space. We wish to associate quantum states with 
our classical initial conditions. The "most classical" state of a quantum system is a 
coherent state, so it would be natural to associate a point in our four-dimensional 
classical phase space with a product of spin coherent states. These states, however, 
have support on several subspaces with different values of F z , and thus correspond 
to a distribution of classical surfaces of sections. To avoid this complication, we 
project our coherent states into the Mp = subspace, and then renormalize them, 
as described in Eq. (3.13). This gives us a pure state, which though no longer 
separable, typically has a low entanglement and is localized around a point in the 
classical phase space in the relevant difference angles. 
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Figure 3.6: Dynamically generated entanglement as a function of the number of 
applications of the Floquet map. (a) Mixed phase space (a = 3/2,(3 = vr/2), regular 
initial condition: \ip Q ) = \I,I)\J,—J). (b) Mixed phase space (a = 3/2,(3 = n/2), 
regular initial condition: \i(j ) = \S9 = n/10, 5(f) = 537r/30). (c) Mixed phase space 
(a = 3/2,(3 = 7r/2), chaotic initial condition: \ip ) = \I , —I)\J, J) (d) Globally 
chaotic phase space (a = 6, (3 = vr/2), chaotic initial condition: |^o) = \S6 = 
tt/2,5<(> = vr/3). The solid line gives the value expected from random states in the 
Hilbert space, Eq. (3.28). 
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The time-evolved state after n applications of the Floquet operator to the pro- 
jected coherent state is 



expanded in the Floquet eigenstates, where a& = (k\56,5<f)) is the initial spectral 
decomposition. The Schmidt coefficients are the expansion of this state in the angular 
momentum product basis (the Schmidt basis) giving, 



according to Eqs. (3.7, 3.20). The Shannon entropy of these coefficients gives the 
dynamically evolved entanglement. Figure 3.6 shows this quantum evolution for 
parameters such that the classical evolution is described by a mixed phase space. 
For a coherent state initial condition chosen in the middle of a regular island (\ipo) = 
\I,I)\J,—J) = \59 = —11, 5(f) = 0)), the entanglement rises slowly and oscillates 
between high and low values. For an initial condition in the chaotic sea (|^>o) — |^ — 
7r/2,50 = 7r/3)), the entanglement rapidly rises and saturates to a near constant 
value, with small fluctuations about the quasi-steady state. 




(3.19) 



k 




2 



(3.20) 
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Figure 3.7: Long-time average entanglement as a function of mean coordinate of 
the initial projected coherent state, (a) Mixed phase space: a = 3/2,(3 = 7r/2 (b) 
Globally chaotic phase space: a = 6, (3 = tt/2) 




Figure 3.8: Side-by-side comparison, showing dynamically generated entanglement 
as superb signature of classical chaos in a mixed phase space (a = 3/2,(3 = tt/2). 
(a) Classical phase space, Poincare section, (b) Long-time average entanglement as 
a function of mean coordinate of the initial projected coherent state 
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In order to better explore how the entanglement evolution saturates to a partic- 
ular value, we average over many time steps to find a long-time average of entan- 
glement. We drop the first three hundred steps in order to remove transient effects 
and insure that the dynamics settle into a quasi-steady state, and then average over 
times steps 300-320. By looking at a plot of this average, we can see how it cor- 
relates with initial conditions, a procedure initially carried out for the kicked top 
Hamiltonian by Wang et al. [70]. Figure 3.8 shows remarkably strong correlation 
between structures in the classical mixed phase space and the long-time entangle- 
ment average plot. Chaotic initial conditions generally go to a higher average value 
than regular initial conditions, with the smallest values of entanglement generation 
near the classical fixed points. Additionally, all initial conditions in the chaotic sea 
saturate to nearly the same average entanglement. 

For parameters corresponding to global chaos, we can see that the surface plot is 
very flat (see Fig. 3.8a), with all initial conditions converging to nearly the same long- 
time entanglement average. For the parameters at hand, averaging over all initial 
conditions, the dynamically generated entanglement is Edynam = 5.28, as compared 
to the value E eigens = 4.97 found for the average entanglement of the eigenstates of 
the Floquet map. For the mixed phase space, the value of long-time entanglement 
is flat for initial conditions that correlate with the classical chaotic sea. To find the 
entanglement characteristic of the chaotic initial conditions, we take a grid of coherent 
states across the phase space. Each point in the grid is determined as "regular" or 
"chaotic" by the local Lyapunov exponent of the classical dynamics. For those states 
with positive Lyapunov exponent we evolve according to the Floquet operator and 
calculate the long-time entanglement average, as described above. Weighting these 
values according to the measure on phase space gives us an average entanglement 
of Edynam = 5.08 in the chaotic sea, significantly lower than that for the globally 
chaotic phase space. Below, we interpret these results with statistics of random 
states in Hilbert space and their connection to quantum chaos. 
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3.3.3 Entanglement and random states in Hilbert space 

The numerical studies in Sec. 3.3.2 reveal some empirical facts. When the Floquet 
map corresponds to a fully chaotic phase space, the entanglement of the eigenstates 
are all nearly equal, with an average value independent of the coupling strength and 
below the maximum possible entanglement for the bipartite system. Moreover, the 
dynamically generated entanglement when starting from a projected coherent state 
localized in a chaotic sea saturates to a nearly constant value after a few applications 
of the Floquet map. In a mixed phase space, the amount of entanglement increases 
as the size of the chaotic sea increases. For a completely chaotic space, the value no 
longer changes with coupling strength. This saturation value is different from the 
entanglement seen in the eigenstates. These facts leads us to conclude that the value 
of entanglement generation for chaotic maps is statistical in nature, as emphasized 
by Bandyopadhyay and Lakshminarayan [63, 58], and Scott and Caves [21]. The 
predicted values follow from the theory of random matrices and random states in 
Hilbert space, which we briefly review. 

The random matrix conjecture of quantum chaos states that when the Floquet 
map (in a periodically driven system) classically generates global chaos, the quantum 
operators have many of the statistical properties of a random matrix drawn from an 
appropriate ensemble depending on fundamental symmetries [9] . Systems with time- 
reversal symmetry (and no Kramer's degeneracy), have Floquet maps with many of 
the statistical properties of random matrices chosen from the Circular Orthogonal 
Ensemble (COE) - the set of symmetric unitary matrices with a probability distribu- 
tion defined by the orthogonally invariant Haar measure [78] . Without time- reversal 
symmetry, the Floquet maps have many of the statistical properties of random ma- 
trices chosen from the Circular Unitary Ensemble (CUE) - the set of general unitary 
matrices with a probability distribution defined by the unitarily invariant Haar mea- 
sure [78]. When expressed in the basis of their eigenvectors, both such matrices have 
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the form U = ^ fc expz0fc|fc)(fc|, where the phases are randomly distributed from 
to 2tt with a uniform probability distribution. In the case of COE, the eigenvectors 
are invariant with respect to an anti-unitary operator T; for CUE the eigenvectors 
have no time-reversal invariance. 

We seek to predict the entanglement of state vectors based on statistical argu- 
ments. We can do this by averaging over an appropriate distribution of random states 
in Hilbert space [' 7]. To construct the probability measure for sampling random 
states, we employ a parameterization equivalent to the Hurwitz parameterization of 
random unitaries [79]. Such a measure can be constructed by connecting the vector 
space with a manifold upon which there is a known geometric measure. A normalized 
state in a <i-dimensional complex Hilbert space can be visualized as a point on the 
surface of a hypersphere in a 2d dimensional real space, where for each of the d basis 
vectors in Hilbert space we assign a pair of orthogonal directions that project out 
the real and imaginary parts of the state's probability amplitude. The surface area 
of a differential patch on a hypersphere is then the probability measure for picking 
uniformly distributed random states. The coordinates of a state, parameterized by 
angles on the hypersphere, and the corresponding measure over the space are 



C X , r = COS 01, 




(3.21a) 


ci t i = sin 9i cos 


02, 


(3.21b) 


c n:r = sin X . . . 


sin 2n -2 cos 2n -x, 


(3.21c) 


c n>i = sin 9i . . . 


sin 2n -x cos 2n , 


(3.21d) 


Cd, r — sin 6\ . . . 


Sin 02d-2COS0 2 d-l, 


(3.21e) 


c d>i = sin X . . . 


sin 2d -x, 


(3.21f) 


dX = Nsm 2d ~ 2 


Ox sin 2d - 3 2 ... , 




sin 2 d-2d0id0 2 


• • • d02d-i, 


(3.21g) 



where c„ jr , c n> j are the real and imaginary expansion coefficients in the n basis 
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state, dX is the surface element, and N is a normalization constant. The angles all 
range from (0, 7r) except for the last angle which varies from (0, 2tt). This defines the 
measure for random vectors over the field of complex numbers. 

For random states in a real vector space, the probability measure is the area 
element on a d dimensional hypersphere, with each direction corresponding to a 
basis vector. In this case the coordinates of the state and measure over the space are 



ci <r = cos 9i, (3.22a) 

c n , r = sin 9i . . . sin(6' ri _i) cos(6 n ), (3.22b) 

c d>r — sin#i . . . sin(0<2-i)> (3.22c) 
dX = Nsm d - 2 6 1 sm d - 3 6 2 ..., 

sm6 d _ 2 d9 1 d6 2 ...d6 2d _ 1 . (3.22d) 



This defines the measure for random vectors over the field of real numbers. 

With these measures in hand, we can calculate expected values of entanglement 
of random states in an appropriate ensemble and compare them to the numerically 
predicted results. For large d-dimensional spaces, the variance scales as \j\fd [ ], so 
when the states in question are well-described by the statistics above, we anticipate 
the expectation value to give good predictive power. A well known example is the 
entanglement of a "typical state" picked at random from a d\ £g> d 2 tensor product 
Hilbert space, with no other restrictions of symmetry. The Haar measure average of 
the entanglement over the whole space gives [80, 21, 81] 

d\d,2 \ (i \ 

E dl ®d2= % 2d — ,d<2 - dl ' (3.23) 

k=di+l 2 

For large dimensions, E dl ®d 2 ~ log^i — di/(2d 2 ), which is close to the maximum 
possible value of entanglement, but saturates slightly below. Typical pure states in 
an unconstrained bipartite Hilbert space are highly entangled [81]. For the case at 
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hand, symmetries constrain the accessible Hilbert space. We thus turn to study the 
typical entanglement expected under these conditions. 

Typical entanglement in a subspace 

Our system has an additional symmetry, its rotational invariance around the z-axis. 
This restricts our system so that eigenstates and dynamics take place in subspaces 
with fixed values of F z . Calculation of entanglement within a subspace is generally 
a nontrivial task as there is no simple expression for the entanglement in terms of 
variables that we can average over the Haar measure [82]. In our case, there is a 
happy accident - the uncoupled basis of angular momentum, | J, mj) ® \I, Mp — m j), 
is also the Schmidt basis for all states in the subspace. This implies that we can take 
the fixed Schmidt vectors as the directions that define the space on a hypersphere, 
and thereby employ the same parameterization of the Haar measures as in Eqs. 
(3.21,3.22), where now d is the dimension of the subspace. Note, this would not 
in general be possible for an arbitrary subspace because the entanglement is not 
a simple function of the expansion coefficients in a fixed basis. The key question 
we must address is whether, with respect to the Schmidt basis, the state vector is 
random over the field of real or complex numbers, since the statistical properties of 
these two vector spaces differ, as discussed by Wootters [77]. Once that question is 
answered, one can predict the entanglement based on the expected entropy in the 
Schmidt basis. 

For a state in a fixed F z subspace, expanded in the uncoupled basis, \^/) = 
J2c mj \mj) \ — mj), the entanglement is 



Note that, the Schmidt basis is T-invariant according to time- reversal operator Eq. 
(3.9). Thus, any other vector which is T-invariant will have real expansion coefficients 




(3.24) 



65 



Chapter 3. Dynamically Generated Entanglement and Chaos 



c mj . If the vector is random with respect to this basis, then these real coefficients 
are distributed on the hypersphere according to Eq. (3.22). The contribution of each 
term in the expression for the entanglement given above should be equal, so we can 
shortcut by integrating only the first term, and multiplying by the number of terms, 
d. We normalize by an integral over the measure for that variable. The result for 
T-invariant vectors is 

-J|cos^| 2 log(|cos^i| 2 ) sin d - 2 M#i 
ER ~ d fsin d - 2 M0i 

= H d/2 + log4-2, (3.25) 

where 

H D = 1 + 1/2 + 1/3 + ■■■ + 1/D (3.26) 

is the harmonic series. 

When the state is not T-invariant, its expansion coefficients in the Schmidt basis 
will be complex. For random complex states states, it is useful to first simplify our 
parameterization by specifying the magnitudes of the expansion coefficients in terms 
of the angles on the hypersphere, rather than the real and imaginary parts of the 
expansion coefficients. Our new parameterization and the associated surface element 
are as follows: 

|ci| = cos 9 1 (3.27a) 
\c m \ = sin 6*i . . . sin# m _i cos^m (3.27b) 
\cd\ = sin d\ . . . sin 9d-i (3.27c) 
dX = Nsm 2d - 3 6 1 sm 2d - 5 6 2 ... 

sin cos 6*i . . . cos9d-id9 1 d9 2 ...d92d-i, (3.27d) 

where 9 m now ranges from (0, 7r/2). Since the entanglement for a state in the sub- 
space depends only on the the magnitudes {|c m |}, Eq. (3.24) can be expressed in 
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terms of this parameterization of the manifold. Performing the average, the typical 
entanglement for a complex state, restricted to a F z subspace, is 

- / | cos 0! | 2 log ( | cos 6 1 1 2 ) sin M - 3 Q x cos Q x dQ x 
E c = d 



/ sin 2d - 3 0! cos M0i 
= H d -1. (3.28) 

These averages hold regardless of dimension of the space, though the variance of the 
distribution rapidly narrows as d increases. 

In the limit of large dimensional spaces, we recover the results of Wootters [' 
and Zyczkowski [83] for the entropy of a random state in a real or complex vector 
space, 

E K -»■ log d - 2 + log 2 + 7, (3.29a) 
E c -»■ logd- 1 + 7, (3.29b) 

where 7 ~ 0.577 is Euler's constant. Whereas these expressions give the entan- 
glement of our state in the F z subspace, in general the entropy of a random state 
with respect to a fixed basis is not equal to its entanglement. For example, for 
the full tensor product space, for large dimensional Hilbert spaces with d\ = c?2, 
Edxm-z ~~ >" l°g^i — 1/2) which differs from the Wooters/ Zyczkowski entropy, Eq. 
(3.29), taking d = d\. 

As an aside, we can repeat our calculations for the linear entropy, an entanglement 
monotone. The linear entropy is determined by the purity of the reduced density 
operator of one subsystem, 

S L (p) = 1 - Tr( Pr 2 ed ) = 1 - Am = 1 " E ' c -| 4 ( 3 - 3 °) 

m m 

2 

where A m = \c m \ are the Schmidt coefficients for a state in the subspace. We repeat 
our integrals over the appropriate manifolds and find 

S LR = 1--^-, 5 LC = 1-^— , (3.31) 
rf + 2' ,c d + V K J 
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the same results found by Brown and Viola by different methods [84]. 

Typical entanglement prediction for the kicked coupled-tops 

With the results of Sec. 3.3.3 in hand, we can compare the predictions of the typical 
entanglement of random states to the entanglement found numerically in Sec. 3.3.2. 
Since the system is time reversal invariant without Kramer's degeneracy as shown 
in Sec. 3.2.2, under the random matrix conjecture of quantum chaos, we expect 
the eigenstates of the Floquet operator for globally chaotic classical dynamics to 
be random real states [85]. The eigenstates are restricted to a subspace with fixed 
value of F z , so Eq. (3.25) applies. We consider the F z = subspace with dimension 
d — 2 J + 1. For spin J = 150, one finds = 4.98, in excellent agreement with the 
mean entanglement of the eigenstates for the globally chaotic case, -E e i g ens = 4.97. 

Next we consider the dynamically generated entanglement, starting from a spin 
coherent product state projected into the F z = subspace. The key conjecture, seen 
numerically in prior studies, is that chaotic maps acting on a fiducial state generate 
states with the statistics of random states in Hilbert space, chosen according to the 
appropriate ensemble. However, contrary to prior claims [63], though the Floquet 
operator is a member of the COE, the dynamically generated state is not a random 
real vector in the Schmidt basis. To see this, first note that since the Floquet 
operator is a member of the COE, we know the eigenstates are time- reversal invariant, 
T\k) = \k). However, according to Eq. (3.9), time reversal acting on the dynamically 
evolved state gives 



Thus, the dynamically evolved state is not an eigenstate of the time reversal operator. 
This is true even when the initial state itself is a time-reversal eigenstate (e.g., the 
coherent state at the pole), in which case T\i/j n ) = \if)- n ). 




(3.32) 
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To further put a point on this, consider the state expanded in the Schmidt basis. 
Simplifying our notation, let \m) = \I,mj = m)\J,mj = —m) be a Schmidt vector. 
After n applications of the Floquet operator, \ip n ) = Ylm c "^\ m )- ^ n ^ ne transforma- 
tion from the initial to final vector in this basis, c>m = E ra '^mm' c m'! the matrix 
^■rnm' = e m< ^* (m \ k) (k \m!) is not an orthogonal matrix. It is a random unitary 
matrix with complex entries with respect to the basis of interest - the Schmidt basis. 
The vector Sm is thus a random vector in complex vector space. 

Given the observations above, we expect the dynamically generated entanglement 
to be predicted by the statistics of random complex vectors. This is indeed born out 
in the numerics. For the globally chaotic map, we evolve and average to find the 
quasi-steady state value, as discussed in Sec. 3.3.2. The long-time entanglement 
average is almost independent of the initial coherent state, projected in the F z = 
subspace. For these initial condition Eq. (3.28) predicts E c = 5.28 in good 
agreement with the long-time average value of 5.28. 

In the case of a mixed phase space, we saw that the long-time entanglement 
average was almost constant for initial states localized in the chaotic sea. Clearly, 
this value of entanglement is a statistical property of Hilbert space. Just as the 
quantum dynamics lead to a random state in the entire F z = subspace when the 
classical dynamics are globally chaotic, for a classically mixed phase space, based on 
Percival's conjecture, the quantum dynamics generate a random state in the chaotic 
subspace. The structure of the chaotic sea cannot be described by a simple symmetry, 
so we cannot determine the entanglement of a typical state analytically. However, 
we can filter the eigenstates to determine which are in the chaotic subspace, as 
discussed in Sec. 3.2.2, and sample randomly from a unitarily invariant measure 
over this subspace in order to find the typical entanglement value. In this case 
there is no simple expression for the entanglement as a function of the states, so 
we cannot analytically take the average over the appropriate measure as before. 
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Instead, we generate a large number of random states in the chaotic subspace, and 
find their entanglements. We do this by picking the real and imaginary parts of the 
expansion coefficients with respect to the chaotic eigenstates according to a Gaussian 
distribution. After normalizing, the entanglement is calculated for this state, and 
the process is repeated 100 times. The results are averaged to find an estimate of the 
average entanglement of a random state in the chaotic sea. We find that the average 
entanglement of a random state in the chaotic subspace picked according to the 
measure for complex random vectors is 5.13, in good agreement with the numerically 
determined value of Edynam = 5.08 found in Sec. 3.3.2. Part of this discrepancy is 
likely due to the greater degree of variation of entanglement across the chaotic sea 
in the mixed phase space compared to the relatively flat completely chaotic phase 
space. In addition our filter for determining the members of the chaotic subspace was 
somewhat crude with an ambiguous "grey zone" . We would expect this to improve 
deeper in the semiclassical regime, where Percival's conjecture applies better. 



3.4 Discussion and Summary 



Classical chaotic dynamics lead to ergodic mixing in phase space. Quantum analogs 
of ergodicity have long been considered, including ergodicity of eigenf unctions [86] , 
"spectral chaos" [87] , and increase in entropy associated with the wave function when 
expanded in a fixed (non-stationary-state) basis [88]. Recent numerical studies in- 
dicate that quantum dynamics generated by nonintegrable Hamiltonians generates 
pseudo-random states in a Hilbert space [21, 89]. In that sense, quantum chaotic 
dynamics is the classical analog of ergodic mixing in quasi-steady state, for times 
sufficiently long compared to the transient behavior, but short compared to the 
Heisenberg time or the time when correlations in the pseudo-random matrix appear 
90]. Such a result is not new, having its roots in the random matrix theory con- 
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jecture of quantum chaos [ ] - the typical Hamiltonian of a nonintegrable system 
has the statistical properties of random matrices of an ensemble picked according 
to the symmetries of the system under time reversal. The classic works on the sub- 
ject, however, focus on the properties of the stationary states and spectra - Berry's 
"quantum chaology" [91]. 

The dynamical generation of random quantum states has implications for the 
dynamical generation of entanglement. It is well known that for large dimensional 
bipartite Hilbert spaces, a random state is highly entangled with almost the max- 
imum entanglement allowed by the dimension [81]. As the large dimensional limit 
is equivalent to the h — > semiclassical limit, and to the degree that the quantum 
analogs of chaotic Hamiltonians generate random states, one expects near maximal 
dynamical generation of entanglement in quantum chaos, to a value that is predicted 
by the statistics at hand. This is not to say that regular dynamics (quantum analogs 
of integrable motion) cannot lead to highly entangled states. Indeed, such behavior 
is seen, and has been previously noted in [62]. Regular dynamics, however, shows 
oscillatory behavior, including in the generation of entanglement. Chaotic dynamics, 
by contrast, lead to quasi-steady state behavior, and typically lead to higher values of 
time-averaged entanglement than regular motion. Taken together, these facts imply 
that the long-time average entanglement in a bipartite system should be a strong 
signature of classical chaos, closely associated with ergodicity in the two dynamical 
descriptions. An initial coherent state in the chaotic sea has a support over a large 
number of Hamiltonian eigenstates. This is in contrast to an initial coherent state in 
the regular islands that have a support over relatively fewer number of eigenstates. 
This causes the dynamically generated entanglement, when starting from a projected 
coherent state localized in a chaotic sea, to saturate to a nearly constant value after 
a few applications of the Floquet map [69]. The entanglement dynamics observed 
are thus closely tied to the support of the initial coherent state in the basis of the 
Hamiltonian eigenstates. The entanglement for a coherent state initially placed in 



71 



Chapter 3. Dynamically Generated Entanglement and Chaos 

a regular island shows a quasiperiodic behaviour due to a support over a relatively 
fewer eigenstates. 

Since entanglement is a global property of the total system, it is critical to study 
the chaos in the joint system dynamics rather than chaos in the separate degrees of 
freedom. It is the joint-system dynamics that mixes the two subsystems and leads 
to random states of the bipartite system with statistically predictable entanglement. 
This perspective helps us to understand some previous results, which though pre- 
dicted analytically and/or numerically, appear to be paradoxical or raise questions 
about the connection between chaos and the dynamical generation of entanglement. 
For example, the results of [61] show that in the case of coupled kicked-tops, when 
the individual subsystems are strongly chaotic but weakly coupled, the rate of gen- 
eration of entanglement decreases with increased chaoticity. This can be understood 
by the fact that in this regime the chaotic mixing between subsystems is suppressed 
due to the increasing mismatch between the time scale governing individual top dy- 
namics and the time scale governing coupling between them. Indeed, Tanaka et 
al. [61] explained this in terms of rapid "dynamical averaging" that washes out the 
correlations that determine the rate of entanglement generation. 

In another example mentioned in the introduction, Demkowicz-Dobrzaski and 
Kus noted that in a highly-regular regime of kicked tops (k = 0.01 chaoticity param- 
eter in the standard notation), the rate of entanglement generation was anomalously 
large. This result, however, is understood by noting that in addition to weak non- 
linearity k for individual tops, the coupling between tops was weak and equal to 
the nonlinear strength, e = 0.01. In that case, all time scales in Eq. 3.1 are of the 
same order, and given the lack of integrability of the total Hamiltonian, H totSu \, we 
expect the global dynamics to be highly chaotic. Indeed, the fact the entanglement 
saturated to the same value seen for chaotic tops without oscillation indicates that 
we reach the entanglement level of a random state in the joint Hilbert space. Our 
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central conclusion is thus that chaos in the subsystem is not a strong indicator of the 
dynamical generation of entanglement, but rather chaos in the joint dynamics of the 
coupled degrees of freedom. The amount of entanglement generated in quasi-steady 
state is statistically predicted by the typical entanglement of a random state in the 
chaotic sea. 

We have studied the relationship between entanglement and chaos for a system 
of isotropically coupled tops in which one of the tops receives a periodic kick around 
a fixed axis. Here the chaos and entanglement arise from the same coupling mech- 
anism removing any ambiguity between chaos in the subsystem vs. chaos in the 
joint-system dynamics. The results reported here give further evidence of the fact 
that chaotic systems take quantum initial conditions to pseudo-random quantum 
states, and that the high long-time entanglement average of states undergoing quan- 
tum chaotic dynamics is just that of a typical state in the Hilbert space. We see the 
confirmation of this picture in the excellent agreement between the properties of en- 
sembles of quantum states and the numerical results for the eigenvector statistics and 
long-time entanglement average for the completely chaotic system. This approach 
was also found to be highly flexible, applying to subspaces and mixed phase spaces. 
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Chapter 4 



Quantum Signatures of Chaos in 
Quantum Tomography 



At a fundamental level chaos represents unpredictability, so this seems at odds with 
the goal of gaining information in order to estimate an unknown quantum state. 
However, on the flip side, this unpredictability represents the potential information 
to be gained in an estimation process. If everything is predicted and known, we learn 
nothing new. The missing information in deterministic chaos is the initial condition. 
A time history of a trajectory at discrete times is an archive of information about 
the initial conditions given perfect knowledge about the dynamics. Moreover, if the 
dynamics is chaotic the rate at which we learn information increases due to the 
rapid Lyapunov divergence of distinguishable trajectories and we expect unbiased 
information because of the ergodic mixing of phase space. That is, if the information 
is generated by chaotic dynamics, the trajectory is random, and all initial conditions 
are equally likely until we invert the data and discover the initial state. 

As in the case of deterministic chaos, the missing information in quantum to- 
mography is the initial condition. We try to accurately model all of the quantum 
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dynamics occurring in the system, and then use the measurement time history to 
give us information about the initial quantum state. The dynamics is informationally 
complete if the time history contains information about an arbitrary initial condi- 
tion. Our goal is to characterize and quantify the performance of tomography, when 
the dynamics driving the system is chaotic in the classical limit. We use this to 
draw comparisons between the role played by regular and chaotic dynamics in the 
information gain in the tomography procedure. 

The remainder of this chapter is organized as follows. In Sec. 4.1, we review the 
protocol for quantum tomography via continuous measurement [27]. We also review 
the essential features of the "quantum chaology" associated with the kicked top [9], 
which is the dynamics employed to drive the system for the purpose of tomography. 
Section 4.2, the heart of this chapter, studies the performance of tomography as the 
system dynamics becomes increasingly chaotic. We perform numerical calculations 
of the average fidelity obtained for an ensemble of random states as a function of 
the chaoticity of the driving dynamics. We then explain these results in terms of 
the information gain in tomography. In section 4.3, we make predictions for the 
information gain using random matrix theory in the fully chaotic regime and show 
a strong agreement between the two. Our results are discussed and summarized in 
Sec. 4.4. 



4.1 Background 

4.1.1 Tomography via weak continuous measurement 

The protocol for quantum tomography via continuous measurement [27, 94] is as 
follows: We are given an ensemble of N, noninteracting, simultaneously prepared 
quantum systems in an identical but unknown state given by the density matrix 
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Pq. Our goal is to determine po by continuously measuring an observable Oq. The 
ensemble is collectively controlled and probed in a time-dependent manner to obtain 
an "informationally complete" continuous measurement record. In order to achieve 
informational completeness, when viewing in the Heisenberg picture, the set of mea- 
sured observables should span an operator basis for p Q . For a Hilbert space of finite 
dimension d, and fixing the normalization of po, the set of Hermitian operators must 
form a basis of su(d). The measurement record is inverted to get an estimate of 
the unknown state. Laboratory Realization of such a record is intimately tied to 
controllability, i.e., designing the system evolution is such a way as to generate ar- 
bitrary unitary maps. While it is desirable to obtain an informationally complete 
measurement record, we shall see that we can obtain high fidelity in tomography in 
some cases even when this is not the case [92]. 

In an idealized form, the probe performs a QND measurement that couples uni- 
formly to the collective variable across the ensemble and measures Y^n=j ■ For a 
strong QND measurement, quantum backaction will result in substantial entangle- 
ment between the particles. For a sufficiently weak measurement, the noise on the 
detector (shot noise of a laser probe) dominates the quantum fluctuation intrinsic 
to the measurement outcomes of the state (projection noise). In this case, we can 
neglect the backaction on the quantum state, and the ensemble remains factorized. 
Then we can write the measurement record obtained as, 

M=Tr(O p(t)) + aW(t), (4.1) 

amplified by the total number of atoms. Here aW(t) is a Gaussian-random variable 
with zero mean and variance a 2 , which accounts for the noise on the detector. 

In order to obtain a measurement record that can be inverted to reconstruct an 
estimate of the initial state, one must drive the system by a carefully designed dynam- 
ical evolution that continually maps new information onto the measured observable. 
In order to do so, the system is manipulated by external fields. The Hamiltonian 
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of the system, H(t) = H[4>i(t)], is a functional of a set of time dependent control 
functions, <j>i(t), so that the dynamics produces an informationally complete measure- 
ment record Ai. Since our goal is to estimate the initial state from the measurement 
record and the system dynamics, we will work in the Heisenberg picture. Rewriting 
eq. (4.1) in the Heisenberg picture, we get 

M = Tr(O(t)p ) + <rWi. (4.2) 

We sample the measurement record at discreet times so that 

Mi = Tr(Oip ) + &Wi. (4.3) 

Thus, the problem of state estimation is reduced to a linear stochastic estimation 
problem. 

The goal is to determine po given {.M;} for a well chosen {Oi} in the presence of 
noise {Wi}. We use a simple linear parametrization of the density matrix 
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where d is the dimension of the Hilbert space, r a are d — 1 real numbers (the 
components of a generalized Bloch vector), and {E a } is an orthonormal Hermitian 
basis of traceless operators. We can then write Eq. 4.3 as 

d' 2 -i 

M,= J2 raTriOiEa) + aW l} (4.5) 

a=l 

or, in the matrix form as, 

M = Or + trW, (4.6) 

which in general is an overdetermined set of linear equations with d 2 — 1 unknowns 
r = (r 1 ,...,r (i 2_ 1 ). 
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The conditional probablity distribution for the random variable M, given the 
state r, is the Gaussian distribution 

P(M|r)ocexp( ^(M — Or) T (M — Or)) (4.7) 

We can use the fact that the argument of the exponent in Eq. 4.7 is a quadratic 
function of r to write the likelihood function (ignoring any priors) 

V(r\M) cc exp(-^(r - r A/L ) T (r - r ML )) (4.8) 

is a Gaussian function over the possible states r centered around the most likely 
state, Yml-i with the covariance matrix given by C = a 2 {0 T 0)~ 1 . The uncontrained 
maximum liklihood solution is given by 

r ML = (d T 6)- 1 6 T M. (4.9) 

The measurement record is informationally complete when the covariance matrix has 
full rank, d 2 — 1. If the measurement record is incomplete and the covariance matrix 
is not full rank, we replace the inverse in Eq. 4.9 with the Moore- Penrose pseudo 
inverse [93]. The eigenvectors of C -1 represent the orthogonal directions in operator 
space that we have measured up to the final time, and the eigenvalues determine the 
uncertainty, or signal-to-noise associated with those measurement directions. 

When we have an incomplete measurement record, or in the presence of noise, the 
unconstrained maximum likelihood procedure does not give a density matrix that 
corresponds to a physical state. The estimated density matrix might have negative 
eigenvalues. We fix this by finding a valid density matrix that is "closest" to Pml, 
the density matrix obtained by the unconstrained maximum likelihood procedure. 
We use the inverse of the covariance matrix, C -1 , as the cost function or metric to 
measure the distance between p ML and p, the new estimate. The metric is defined 
as 

II r ML - f || 2 = (r ML - v)C- l {v ML ~ f). (4.10) 
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This metric can be justified in the following way. The inverse of the covariance ma- 
trix, C _1 , encodes the nature of information obtained by our procedure. A small 
eigenvalue of C" 1 , implies that we have a low signal to noise ratio associated with the 
measurement of the eigen-operator or more uncertainty in that particular direction 
in the operator space. By defining the cost function as in Eq. 4.10, we take into 
account the fact that different directions in operator space have been measured un- 
equally. The metric in Eq. 4.10 makes sure that during the numerical optimization 
while finding a positive matrix, p, the uncertain components can be adjusted with 
more freedom that the more certain ones thereby maintaining faithfulness with the 
measurement record. 

Thus we see, the key quantity that quantifies the information gain about the state 
is the inverse of the covariance matrix, C -1 the multivariate Gaussian associated with 
the measurement record. Using C _1 , we will be constructing metrics for information 
gain in tomography and see how these metrics behave as we change the control 
dynamics of our experiment from regular to chaotic. 

4.1.2 The kicked top 

How does the presence of chaos in the control dynamics influence our ability to 
perform tomography? In order to address this question, we chose the "kicked top" 
dynamics [ I] as the paradigm to explore quantum chaos in tomography. The Hamil- 
tonian for the kicked top (after setting h — 1) is given by 



Here, the operators, J x , J y and J z are the angular momentum operators obeying the 
commutation relation [J i: Jj] = ie^Jk The first term in the Hamiltonian describes 




n=— oo 



(4.11) 
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a precession around the x axis with an angular frequency - and the second term 
described a periodic sequence of kicks separated by time period r. Each kick is an 
impulsive rotation about the z axis by an amount proportional to J z . Choosing the 
external field to act in delta kicks allows us to express the Floquet map (transfor- 
mation after one period) in a simple form of sequential rotations, 

U T = e = ^e~ iaJ ' 1 (4.12) 

where a and A are related to p and k in terms of the kicking period. The classical 
map can be obtained by considering the Heisenberg evolution of the expectation 
values of the angular momentum operators in a familiar way [9]. 

The classical dynamics consists of the motion of a unit spin vector on the surface 
of the sphere. The z-component of a spin and the angle 0, denoting its orientation in 
the x-y plane, are canonically conjugate, and thus the spin constitutes one canonical 
degree of freedom. The classical dynamical map has the same physical action as 
described above in the quantum context - precession of the spin around the x axis 
with an angular frequency a followed by an impulsive rotation around the z axis by 
an amount proportional to J z and a proportionality constant A. In our analysis, we 
fix a = 1.4 and choose A to be our chaoticity parameter. As we vary A from to 
7, the dynamics change from highly regular to completely chaotic. Since the total 
magnitude of the spin is a constant of motion, our classical map is two dimensional. 
We visualize the phase by plotting the z and y components of motion after every 
application of the dynamical map. 

In order to explore the connection between the degree of chaos and the informa- 
tion gain in tomography, we will consider the kicked top dynamics in four different 
regimes, with varying degrees of chaos. Figure 4.1 shows four different regimes of 
classical dynamics. With the parameters a = 1.4, A = 0.5 (Fig. 4.1a), the dynamics 
are highly regular. When a = 1.4, A = 2.5 (Fig. 4.1b), we see a mixed space with 
chaotic and regular regions of comparable size. The parameters a = 1.4, A = 3.0 (Fig. 
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Figure 4.1: Phase space plots for the kicked top in four regimes, (a) Regular phase 
space: a = 1.4, A = 0.5, (b) Mixed phase space: a = 1.4, A = 2.5, (c) Mostly chaotic: 
a = 1.4, A = 3.0, (d) Fully chaotic phase space: a = 1.4, A = 7.0. The figures depict 
trajectories on the southern hemisphere (x < 0) of the unit sphere where X = y, 

Y = ^ and Z = ^? and we take the limit j — > oo to get the classical limit as in [9]. 



4.1c), give a phase space that has mostly chaotic regions and finally, a = 1.4, A = 7.0 
gives a completely chaotic phase space (Fig. 4. Id). 
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4.1.3 Metrics to quantify information gain 

Our protocol for quantum tomography via continuous measurement of a driven sys- 
tem [: '] gives us a window into the complexity of quantum dynamics and its rela- 
tionship to chaos. The experimental implementation of tomography by continuous 
measurement provides a nice platform for exploring these ideas in the laboratory [28] . 
Quantum tomography deals with extraction of information about an unknown quan- 
tum state through measurements. In our attempt to study dynamical chaos under 
this paradigm, we define metrics to quantify this information gain. These metrics 
allow us to study the ability of our control dynamics to generate a sufficiently high 
signal-to-noise ratio in different directions of the operator space. As we shall see, 
these metrics elucidate the connection between the degree of chaos and the fidelities 
obtained in tomography. We can quantify the information gain in a number of ways. 

1) If we normalize the eigenvalues of the inverse of the covariance matrix, then as 
a probability distribution, its Shannon entropy, E is a measure of how evenly we have 
sampled all the directions in the operator space. We reach the maximum entropy 
when we have measured all directions in the space of matrices equally, E max = 
\og(d 2 — 1). This is the most unbiased measurement we can implement that will lead 
to the highest fidelities, on average, for a random state. 

2) The Fisher information sets the Cramer-Rao bound on our uncertainty about 
the parameters that define the density operator. In general, for a multivariate pa- 
rameter estimation problem, the Cramer-Rao bound gives 

Cov e (T(X)) > F" 1 (4.13) 

where, matrix inequality A > B is understood to mean that the matrix A — B is 
positive semidefinite. T(X) is the unbiased estimator of the multivariate parameter, 
6 = 62, Od) an d Covg(T(X)) is the covariance matrix of a set of unbiased 
estimators for the parameters 9. It quantifies the error in our estimation process. F 
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is the multivariate generalization of the Fisher information [32] 

F mn = EkJ-log/fo 6)^-\ogf(x; 9)). (4.14) 
ov m ou n 

For a multivariate Gaussian, the Fisher information matrix is equal to the inverse of 
the covariance matrix 

F = CT 1 (4.15) 
and thus the Cramer-Rao bound reads 

Cov e (T(X)) > C (4.16) 
We consider the basis in which F and hence C~ l is diagonal. 

F' = UFU T . (4.17) 

Such a transformation is provided by U composed from the eigenvectors of C. In 
this representation, the estimate of the newly transformed parameters fluctuate in- 
dependently of each other. This suggests the possibility to form a single number 
quantifying the performance of the tomography scheme as a whole by adding those 
independent errors 

e>Tr(C) (4.18) 

Thus, jy^y, which is the collective Fisher information (F. I.), serves as a measure of 
the amount of information about the parameter 9 that is present in the data. 



4.2 Fidelity in quantum tomography as a signa- 
ture of chaos 

How does the presence of chaos in the control dynamics influence our ability to 
perform tomography? In order to address this question, we chose the "kicked top" 



83 



Chapter 4. Quantum Signatures of Chaos in Quantum Tomography 

dynamics [ ] as the paradigm to explore quantum chaos in tomography. While one 
can show that repeated application of such a single unitary, like the kicked top, does 
not generate an informationally complete record [92] , one can still obtain sufficiently 
high fidelities taking into account the positivity constraint in the numerical opti- 
mization. In [92] , it was shown that such a repeated application of the "kicked top" 
dynamics for the purpose of tomography is equivalent to the repeated application of 
a single random unitary as far as the final fidelities of reconstruction are concerned. 
The arguments made in [! !] were more "kinematic" in nature as they mostly dealt 
with the properties of random matrices and how close such a dynamics driven by a 
single random unitary comes to obtaining an informationally complete measurement 
record. In this chapter we explore, even more thoroughly, the role of chaos in infor- 
mation gain in tomography. Our approach is more "dynamical" in the sense that 
we are not only interested in the final fidelities obtained, but also how the rate of 
information gain is influenced by the degree of chaos in the dynamics. 

We are now ready to explore the role of chaos in the performance of tomography. 
Throughout this section, we consider spin J = 10, which is sufficiently large that a 
minimum uncertainty spin coherent state can resolve features in the classical phase 
space. Figure 4.2 shows the average fidelity of reconstruction of 100 states picked at 
random according to the Haar measure as a function of the number of applications of 
the kicked top map, and for different values of the chaoticity parameter. We see that 
the rate of increase in fidelity increases with the degree of chaos. The final fidelity 
achieved after a fixed number of kicks is also correlated with the degree of chaos. 
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Figure 4.2: Fidelity of reconstruction as a function of the number of applications of 
the kicked top map. The fidelity calculated as the average fidelity of reconstruction 
of 100 states picked at random according to the Haar measure. The parameters of 
the kicked top are as described in the text, with a fixed. We show the fidelity for 
different choices of the chaoticity parameter. Both the rate of growth and the final 
value of the fidelities are increased with higher values of A. 
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We can understand the above results by studying the information gain in to- 
mography as a function of the degree of chaos in the control dynamics. Figure 4.3 
shows the behavior of the entropy E of the covariance matrix, as defined above, as a 
function of the number of applications of the kicked top map, and for different values 
of the chaoticity parameter. We see that the rate of entropy increase for short times, 
Fig. 4.3a, is correlated with the degree of chaos present in the control dynamics. The 
asymptotic value of the entropy reached also increases with the chaoticity parame- 
ter. Chaotic dynamics provides a measurement record with a large signal-to-noise 
ratio in all the directions in the operator space. Increase in the chaoticity parameter 
results in an increasingly unbiased measurement process that will yield high fideli- 
ties for estimating random quantum states. Figure 4.3a shows the behavior of the 
entropy at short time scales, while we see asymptotic behavior in Fig. 4.3b. The 
collective Fisher information, , tells us about the amount of information our 
measurement record contains about the parameters that define the density matrix. 
Figure 4.4 shows the behavior of the Fisher information as a function of the num- 
ber of applications of the kicked top map, and for different values of the chaoticity 
parameter. We see that the rate of increase of the Fisher information is correlated 
with the degree of chaos present in the control dynamics. As our dynamics become 
increasingly chaotic, we obtain higher values for the Fisher information at a given 
time. We expect the Fisher information to be correlated with the average fidelities 
of estimation for an ensemble of random states. 
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Figure 4.3: The Shannon entropy of the normalized eigenvalues of the inverse of 
covariance matrix as a function of the number of applications of the kicked top map: 
(a) Short time behavior (b) long time/asymptotic behavior. The parameters are as 
described in the text. 



87 



Chapter 4. Quantum Signatures of Chaos in Quantum Tomography 




Figure 4.4: The Fisher information of the parameter estimation in tomography as a 
function of the number of applications of the kicked top map: (a) Short time behavior 
(b) long time/asymptotic behavior. The parameters are as described in the text. 
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4.3 Signatures of chaos : Information gain in the 
fully chaotic regime and random matrix the- 
ory 

As discussed in the previous chapters, a central result of quantum chaos is the con- 
nection with the theory of random matrices [9]. In the limit of large Hilbert space 
dimensions (small h), for parameters such that the classical description of the dy- 
namics shows global chaos, the eigenstates and eigenvalues of the quantum dynamics 
have the statistical properties of an ensemble of random matrices. The appropri- 
ate ensemble depends on the properties of the quantum system under time-reversal 
symmetry [9]. We thus seek to determine whether there exists an anti-unitary (time 
reversal) operator T that has the following action on the Floquet operator, 



2 



TU T T- 1 = XJ\ = e iaJ *e^f. (4.19) 
Considering the generalized time reversal operation, 

T = e iaJ *K, (4.20) 
where K is the complex conjugation operator. It then follows that 

TU T T- 1 = (e iaJ *K) fe^e-^A (Ke~ iaJ *) (4.21) 



e 



iaJ x 



+i\j 2 



so the dynamics is time- reversal invariant. Moreover, T 2 = 1, so there is no Kramer's 
degeneracy. Given these facts, for parameters in which the classical dynamics is 
globally chaotic, we expect the Floquet operator to have the statistical properties of 
a random matrix chosen from the circular orthogonal ensemble (COE). 
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When the system is driven by dynamics that are completely chaotic, we expect 
the information gain and the fidelity to follow the predictions from random matrix 
theory. Figure 4.5 shows the behavior of the fidelity, Shannon entropy and the Fisher 
information of the inverse of the covariance matrix as a function of the number of 
applications of the kicked top map (the blue line) and compares it with the Shannon 
entropy a typical random unitary picked from the COE (the green line). We see a 
remarkable agreement between our predictions from random matrix theory and the 
entropy calculation for the evolution by a completely chaotic map. 
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Figure 4.5: Comparison between the tomography performed by the repeated appli- 
cation of kicked top in the fully chaotic regime (the blue line) and that by a typical 
random unitary picked from the COE (the green line), (a) The average fidelity of 
reconstruction of 100 states picked at random according to the Haar measure, (b) 
The Shannon entropy of the normalized eigenvalues of the inverse of covariance ma- 
trix as a function of the number of applications of the map. The dotted line gives 
the upper bound on the entropy E max = log(cP — 1). 
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We test our predictions from the random matrix theory for chaotic maps without a 
time reversal symmetry. For example, another type of the "kicked top" map without 
time reversal symmetry is given by 

JJ _ g— iAi — J%—iai J x Q-i>~2- Jy—ia-2. J y ^—1X3 J^—ia.3 J z (a 22) 

In Fig. 4.6, we repeat the above calculations for this map. In this case, the appropri- 
ate random matrix ensemble is the CUE. We see an excellent agreement between the 
behavior of the Shannon entropy as predicted by the random matrix theory and that 
for the evolution by a completely chaotic map without the time reversal symmetry 
[95]. 
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Figure 4.6: Same as Figure 4.5 but for a kicked top without time reversal invariance 
(Eq. 4.22). In this case the results are well predicted by modeling the dynamics by 
random matrices sampled from the CUE. 
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When all the eigenvalues of the inverse of the covariance matrix are equal, we 
have an upper bound on the entropy, E max = \og(d 2 — 1). Figures 4.5 - 4.6 compare 
the entropy values achieved by the repeated application of the same unitary (time 
reversal invariant or otherwise) to E max . We see that we fall significantly short of 
Emax by such a procedure. 

So far, we have considered the application of the same unitary matrix over and 
over again to obtain the measurement record. However, this does alone give us an 
informationally complete measurement record and high fidelities are reached only 
when we make use of the positivity constraint. On the other hand, we can consider 
application of a series of different random unitaries, randomly chosen. In that case 
we expect to rapidly reach an informationally complete set and thus rapidly gain 
information about tomography. In Fig. 4.7, we plot the fidelities, Shannon entropy 
and Fisher information achieved by applying a different random unitary, picked from 
the unitarily invariant Haar measure, and compare it with the results obtained by the 
repeated application of the same unitary (picked from the COE and CUE). We also 
see that we reach the upper bound, E max , asymptotically, by this method. Indeed, 
an application of a different random unitary is the most unbiased estimation we can 
hope to perform. 
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Figure 4.7: Comparison between tomography performed by applying a different ran- 
dom unitary at each time step, picked from the unitarily invariant Haar measure 
(magenta line) and that by a repeated appplication of a random unitary picked from 
the COE (blue line) and the CUE (red line) . The dotted line gives the upper bound 
on the entropy, E max = log(oP - 1) 
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4.4 Discussions and Summary 

Dynamics sensitive to the initial conditions will reveal more information about them 
as one observes the system trajectory in the course of time. Classically chaotic dy- 
namics, being sensitive to the initial conditions, generates this unpredictability or 
information to be gained about the initial coordinates of the trajectory. Similarly, 
we found that the rate at which one obtains information about an initially unknown 
quantum state in quantum tomography is correlated with the extent of chaos in the 
system. This is a new quantum signature of classical chaos that we have found. In 
fact, our results can be regarded as signatures of chaos in closed quantum systems 
undergoing unitary evolution. We have been able to quantify the information gain 
using the Fisher information associated with estimating the parameters of the un- 
known quantum state. When the system is fully chaotic, the rate of information gain 
agrees with the predictions of random matrix theory. 

At its core, our approach is akin to the Kolmogorov- Sinai (KS) entropy measure 
of chaos [. 2]. Incomplete information about the initial condition leads to unpre- 
dictability of a time history. In the presence of classical chaos, in order to predict 
which coarse-grained cell in phase space a trajectory will land at a later time, we 
require an exponentially increasing fine-grained knowledge of the initial condition. 
The KS entropy is the rate of increase, and is related to the positive Lyapunov ex- 
ponents of the system. Is there a meaningful quantum definition of KS entropy? In 
order to predict the measurement record with a fixed uncertainty we learn more and 
more about the initial condition. Is the rate at which we obtain this information 
exponentially fast when the system is quantum chaotic? Does this converge to the 
classical Lyapunov exponents in the limit of large action (small H)7 There are many 
important subtleties in these questions. As we gain more and more information, 
eventually quantum backaction becomes important in the measurement history. The 
number of copies we have and the shot noise on the probe limits the ultimate res- 
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olution with which we can deduce the quantum state. Unlike classical dynamics, 
we can never consider infinite resolution, even in principle. The quantum resolution 
is limited by the size of h. As the dimension of the Hilbert space increases, and 
hence the effective h decreases, we expect to see an even sharper difference in the 
information gain as a function of chaoticity. In the limit when d, the dimension of 
the Hilbert space, becomes infinity, we expect the rate of information gain to be inti- 
mately related to the classical Lyapunov exponents. How does all this translates into 
a quantum definition of KS entropy is an important subject of further investigation. 

Still further, in principle we never have perfect knowledge of the dynamics. An- 
other hallmark of chaos is hypersensitivity to perturbations [ ] . How does this funda- 
mentally limit our ability to perform quantum state reconstruction when the system 
is sufficiently complex, and the equivalent dynamics is chaotic. Though quantum 
systems show no sensitivity to initial conditions, due to unitarity, they do show a 
sensitivity to parameters in the Hamiltonian [96]. This poses interesting questions 
for quantum tomography and, more interestingly, for quantum simulations. Under 
what conditions are the system dynamics sensitive to perturbations and how does 
this effect our ability to perform quantum tomography? Under what conditions does 
the underlying quantum chaos affect our ability to accomplish quantum simulations 
in general? We hope to answer these questions in our future work. 



97 



Chapter 5 



Interpreting Quantum Discord 
through Quantum State Merging 

5.1 Introduction 

The characterization of the resources behind the enhancements and speedups pro- 
vided by quantum mechanics over best known classical procedures is one of the most 
fundamental questions in quantum information science. This has generally been 
expressed in terms of quantum entanglement [ 7]. There are, however, quantum 
processes which provide an exponential advantage in the presence of little or no en- 
tanglement [98, 99]. In the realm of mixed-state quantum computation, quantum 
discord [2, 50] has been proposed as a resource [3] and there has been continued 
progress [100]. It has also been shown to be a resource in quantum state discrimina- 
tion [101, 102] and quantum locking [103]. 

Quantum discord was originally suggested as a measure of quantumness of cor- 
relations [2], and has since been studied in variety of systems and settings [122, 123, 
124]. The initial motivation for its definition arose in the context of pointer states 
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and environment-induced decoherence [125]. It has since been related to quantum 
phase transitions [ 26] and the performance to quantum and classical Maxwell's 
demons [127]. Though satisfactory from a physical perspective, a benchmark for ac- 
cepting some quantity as a resource in quantum information science is that it appears 
as the solution to appropriate asymptotic information processing tasks. It is this in- 
formation theoretic interpretation that has been lacking for quantum discord, and 
which we now provide in this dissertation. This also addresses a more fundamental 
dichotomy in quantum information science, where resources and their manipulations 
can have both a thermodynamic and an information theoretic interpretations inde- 
pendently, which are not intuitively or mathematically reconciled. 



5.2 Quantum State Merging And Discord 

Consider a party Bob having access to some incomplete information Y, and another 
party Alice having the missing part X. We can think of X and Y as random variables. 
If Bob wishes to learn X fully, how much information must Alice send to him? 
Evidently, she can send H(X) bits to satisfy Bob. Here H(.) is the Shannon entropy 
associates with the random variable under consideration. However, Slepian and 
Wolf showed that she can do better, by merely sending H(X\Y) = H(X,Y) — 
H(Y), the conditional information [128]. Since H(X\Y) < H(X), Alice can take 
advantage of correlations between X and Y to reduce the communication cost needed 
to accomplish the given task. 

The quantum state merging protocol is the extension of the classical Slepian- Wolf 
protocol into the quantum domain where Alice and Bob share n copies of a quantum 
state, expressed in the form of a tensor power as with each party having the 
marginal density operators p® n and p® n respectively. Let \^abc) be a purification 
of Pab- Assume, without loss of generality, that Bob holds C. The quantum state 
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merging protocol quantifies the minimum amount of quantum information that Alice 
must send to Bob so that he ends up with a state arbitrarily close to \^)b^bc-> B' 
being a register at Bob's end to store the qubits received from Alice. It was shown 
that in the limit of n — > oo, and asymptotically vanishing errors, the answer is given 
by the quantum conditional entropy [149, 150], S(A\B) = S(A,B) - S(B). When 
S(A\B) is negative, Bob obtains the full state with just local operations and classical 
communication, and distills —S(A\B) ebits with Alice, that can be used to transfer 
additional quantum information in the future. 

It is to be noted that quantum discord can be nonzero even for separable quantum 
states. The state merging protocol does in fact work in such a scenario as was shown 
by Horodecki et. al [114]. In such a case, the entanglement needed for state merging 
is provided by sharing the right amount of EPR pairs between Alice and Bob. The 
cost of sharing these EPR pairs is properly taken into account when one computes 
the net cost of quantum communication needed to accomplish state merging. We 
quote from Horodecki et. al. [114]: 

Theorem 2 (Quantum State Merging). For a state pab shared by Alice and Bob, the 
entanglement cost of merging is equal to the quantum conditional entropy S(A\B) = 
S(AB) — S(B) in the following sense. When the S(A\B) is positive, then merging is 
possible if and only if R > S(A\B) ebits per input copy are provided. When S(A\B) 
is negative, then merging is possible by local operations and classical communication, 
and moreover R < —S(A\B) maximally entangled states are obtained per input copy. 

Although it is true that one needs entanglement for state merging (and thus one 
needs to share extra EPR pairs) in the case when we have a separable state having 
nonzero discord, we demonstrate that even in such a scenario, discord is equal to the 
mark up in the amount of entanglement needed for state merging after one of the 
parties (B) makes a measurement. 

In both scenarios (independent of whether S(A\B) is positive or negative), the 
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cost of state merging is more expensive when one destroys the quantum part of 
correlations (which are quantified by discord) , by making measurements on one sub- 
system B. Thus we see, state merging and hence our interpretation is valid in general 
including the cases when we have a separable state with nonzero discord. 

5.2.1 Quantum operations increase the cost of state merging 

An intuitive argument for our interpretation of quantum discord begins with the 
strong subadditivity theorem [112], which states that [150] 

S(A\B,C) <S(A\B). (5.1) 

From the point of view of the state merging protocol, the above has a very clear 
interpretation: having more prior information makes state merging cheaper. Or in 
other words, throwing away information will make state merging more expensive. 
Thus, if Bob discards system C, it will increase the cost of quantum communication 
needed by Alice in order to merge her state with Bob. Our intent here shall be to 
relate this increase in the cost of state merging to quantum discord between A and 
B. 

To that end, we expand the size of the Hilbert space so that an arbitrary mea- 
surement (with forgetting) can be modeled by coupling to the auxiliary subsystem 
and then discarding it. This permits us to apply strong subadditivity to the problem 
in question. We assume C is initially in a pure state |0), and a unitary interac- 
tion U between B and C. Letting primes denote the state of the system after U 
has acted we have S(A, B) = S(A, BC) as C starts out in a product state with 
AB. We also have I(A : BC) = I (A : B'C). As discarding quantum systems 
cannot increase the mutual information, we get I {A : B') < I (A : B'C). Now 
consider the state merging protocol between A and B in the presence of C. We have 
S(A\B) = S(A) - I (A : B) = S{A) - I (A : BC) = S{A\BC). After the application 
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of the unitary U, but before discarding the subsystem C, the cost of merging is still 
given by S{A'\B'C) = S(A\B). This implies that one can always view the cost of 
merging the state of system A with B, as the cost of merging A with the system BC, 
where C is some ancilla (initially in a pure state) with which B interacts coherently 
through a unitary U. Such a scheme does not change the cost of state merging, as 
shown, but helps us in counting resources. Discarding system C yields 



Computing the mark up in the price in the state merging on discarding information 
gives D = I (A : B) — I {A 1 : B'). This quantity D is equal to quantum discord when 
our quantum operations are quantum measurements maximizing I (A' : B'). Thus, 
discord is the minimum possible increase in the cost of quantum communication 
in performing state merging, with a measurement on the party receiving the final 
state. This also addresses the asymmetry that is inherent is quantum discord. This 
is exhibited operationally in our interpretation since the state merging protocol is 
not invariant under exchanging the parties. 

We now show that the minimum of D over all possible measurements is the 
quantum discord. The state pab, under measurement of subsystem B, changes to 
Pab = J2j PjPA\j® 7I 'j, where {ttj} are orthogonal projectors resulting from a Neumark 
extension of the POVM elements [131]. The unconditioned post-measurement states 
of A and B are 



I {A' : B') < I {A' : B'C) = I (A : BC) = I (A : B) 



(5.2) 



or alternatively, 



S{A'\B') > S{A'\B'C) = S(A\B). 



(5.3) 




j 



3 
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Computing the value of I (A' : B'), we get 



I (A' : B') 



S(A') + S(B') - S(A',B'), 

S(A') + H{p) - {H(p) + 5>£(p^)}, 



3 



S{A)-Y,PiS{PA\ j ). 



(5.4) 



3 



After maximization, it reduces to JJ(pab), as in Eq. (2.28). The reduction to rank 
1 POVMs follows as stated earlier. 

We can also rewrite the expression for D using Eq. (5.3) instead of Eq. (5.2) as the 
increase of the conditional entropy D = S(A'\B') — S(A\B). The above expression 
makes our interpretation even more transparent. Quantum measurements on B 
destroys quantum correlations between A and B. This increases the average cost of 
quantum communication needed by A, to merge her post measurement state with 
B. Since, S(A'\B') = J2jPjS(pA\j) > S(A\B), there is always a mark up in the cost 
of state merging. 

5.3 An Example 

Consider the separable state p AB = (|0) A (0| <g> |0) B (0| + ® l+M+l) / 2, where 

|+) = (|0) + \l))/y/2. The cost incurred by A to merge her state with B is S(A\B) = 
0.399124 ebits, while that after measuring B using the projectors (J ± ax ^ z ) / 2 is 
S(A'\B') = 0.600876 ebits. The markup in the cost of state merging is S(A'\B') — 
S(A\B) = 0.201752 ebits, which is exactly the quantum discord of the state qab- 
Hence, any information lost through the measurement results in making the quantum 
state merging more expensive by exactly the same amount. 
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5.3.1 Properties of Discord 

We can now use our quantum state merging perspective to derive the various proper- 
ties of discord. Since measurements on system B will always result in either discard- 
ing some information or at best preserving the original correlations, we will always 
get a price hike in state merging or at best we can hope to just break even. Hence, 
discord, which is the mark up, will always be greater than or equal to zero [2, 132]. 

Quantum discord of a state is zero if and only if the density matrix is of the 
form pab = YliiViPA\i ® i n the basis which diagonalizes ps- Measuring the 

projectors |Aj)(Aj| and discarding the measurement results on such a state yields 
Pab = ^i-PjPABPj = Pab- Thus, we have a measurement which causes no loss of 
information, and retains all the correlations between A and B. Hence there is no 
mark up in the cost of merging a zero discord state. 

The converse can be seen through the application of strong subadditivity in 
Eq. (6.13). The equality of mutual information, I(A : B), of the initial state and that 
of the state after quantum operations on B, I (A' : B') coincides with the equality 
condition for strong subadditivity. But this is exactly the condition for quantum 
discord [ ] to vanish. Thus a zero mark up in the cost of state merging implies 
zero discord. 

An upper bound on discord is decided by an upper bound on the mark-up we 
can obtain. Since Bob cannot lose more information than there is at his disposal, 
the entropy of the state at Bob's end, S(B), is an upper bound on quantum discord. 

Finally, for pure states, quantum discord reduces to entanglement, and S(A\B) 
= S(A) — I (A : B) = —S(A) < 0. From our perspective, measurement destroys 
all the correlations present between A and B. Though the post-measurement state 
merging of the state of A with that of B occurs at zero cost, they lose the —S(A\B) 
potential Bell pairs that could have been put to some use. This provides us a novel 
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way of measuring entanglement, as the markup in merging a pure state, when B is 
measured. 

In conclusion, this study places quantum discord squarely in the midst of quantum 
informational concepts and opens up the way for its manipulation as a resource 
in quantum information processing. We also hope that our work will serve as a 
stepping stone for a more comprehensive and unified understanding of quantum 
physics, thermodynamics and information theory. 
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Quantum discord as a resource in 
quantum protocols 

Quantum correlations lie at the heart of the mystery of quantum mechanics and 
also serve as a resource for the potential benefits quantum information science might 
provide. However, recent developments suggest that entanglement may not fully 
capture the complete quantum character of a system [2]. Quantum discord aims to 
fill this gap and quantify essentially all the quantum correlations in a quantum state. 

In the previous chapter, we addressed the important question: Does quantum 
discord have a definite physical role in information processing? We provided an affir- 
mative answer by linking quantum discord to the yield of physical tasks in quantum 
information theory. At this point, a natural question arises: Is there a role of quan- 
tum discord in quantum information theory as a whole beyond the state merging 
protocol? This chapter answers precisely this question. 

We show quantum discord to be a resource in quantum information process- 
ing. This is accomplished by proving a relationship between quantum discord and 
the yield of the quantum protocols. Our results are derived by studying the fully 
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quantum Slepian-Wolf protocol[101] - a unification of essentially all bipartite, uni- 
directional and memoryless quantum communication protocols - in the presence of 
environmental decoherence. As examples, we elucidate the significance of quantum 
discord in quantum teleportation, superdense coding, and entanglement distillation. 
Finally, we provide the first quantitative relation between quantum discord and the 
query complexity of quantum computations. 

The key insight to our findings is that quantum measurements and environmental 
decoherence disturb a quantum system in a way that is unique to quantum theory. 
Quantum correlations in a bipartite system are precisely the ones that are destroyed 
by such disturbances, and therefore certain quantum communication protocols be- 
come overloaded by an amount exactly equal to quantum discord. More specifically, 
we showed that discord is the markup in the cost of quantum communication in 
the process of quantum state merging [111], if the system undergoes measurement 
and/or decoherence. At this point, a natural question arises: Is there a role of 
quantum discord in quantum information theory as a whole beyond the lossy state 
merging protocol? We answered this question with a yes. We first observed that 
quantum state merging protocol is a derivative of the more general Fully Quantum 
Slepian Wolf (FQSW) protocol [114, c ] and the closely related "mother" proto- 
col. The mother protocol is essentially a unification of essentially all unidirectional, 
bipartite and memoryless quantum communication protocols like quantum telepor- 
tation, superdense coding and entanglement distillation. A link between discord and 
state merging does indeed suggest a link between discord and the mother protocol 
and hence a possible role of discord in essentially all bipartite, unidirectional and 
memoryless quantum communication protocols. 

This is made possible by comparing the performance of the fully quantum Slepian- 
Wolf (FQSW) protocol [104] in the presence and absence of decoherence and linking 
it to the discord of the state involved. While decoherence is expected to diminish the 
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gain provided by a quantum protocol, we provide, for the first time, a general lower 
bound on the amount of this deterioration. Our bound is only dependent on the 
state involved, independent of the details of the protocol as well as the nature of the 
decoherence. Within the resource framework of quantum Shannon theory [104, 30], 
we couple the performance of the FQSW protocol to the most general environmental 
decoherence to show that quantum discord of the state participating in the protocol 
is the lower bound to the depreciation of the protocol's performance. The FQSW 
protocol - a quantum communication-assisted entanglement distillation protocol - is 
the parent protocol from which all information processing protocols emanate [33]. 
The generality of the FQSW protocol allows us to establish the role of quantum 
discord in the performance of noisy versions of quantum teleportation, super-dense 
coding, and distillation. Finally, interpreting quantum computation as the classical 
capacity of a quantum communication channel [ I I 5] , we also elucidate the role of 
quantum discord in quantum computational processes. 



Our result connecting quantum discord to quantum protocols is subtle as well as 
intriguing. Although it is known that entanglement is often necessary for the success 
of quantum protocols, and that the presence of decoherence affects its performance, 
we have now provided a quantitative result of the amount of such a depreciation. 
We have shown that the amount by which a protocol suffers in the presence of 
decoherence is an inherent property of the quantum states involved. It suggests that 
the choice of the best state for any noisy quantum protocol must be a tradeoff between 
the entanglement and discord of the state involved. Given the non-monotonic relation 
between quantum discord and entanglement in quantum states [110], choosing the 
optimal state for a quantum task is a non-trivial one, though for which our work 
identifies the proper certificate. 
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6.1 The FQSW protocol 

6.1.1 The mother protocol and the quantum information's 
family tree 

So far, we have been discussing quantum communication tasks like teleportation, 
super-dense coding, quantum state merging as independent entities. Is there a uni- 
fied way of looking at various quantum communication tasks? Do various quantum 
communication protocols described above have a common origin. Interestingly, the 
answer is yes. It was shown in [33] that essentially all unidirectional, bipartite and 
memoryless quantum communication protocols are actually siblings originating from 
one "mother". The mother protocol can be seen to provide a hierarchical struc- 
ture to the family of quantum protocols. We will describe the improved version 
of the mother protocol, which is the fully quantum Slepian-Wolf (FQSW) protocol. 
Throughout this chapter, we will use the names "mother" and FQSW to describe 
essentially this same protocol. We have described the FQSW protocol in detail in 
Chapter 2 (section 2.4). In the next section, we examine the FQSW protocol in the 
presence of environmental decoherence. 

6.2 Quantum discord as a measure of coherence 
in the FQSW protocol 

This section contains our main result. The FQSW is essentially a non-dissipative 
protocol in that no information is leaked to the environment in each step of the 
protocol, but any practical implementation of a quantum information protocol will 
be affected by loss and noise. In particular, we will consider loss of information 
and coherence at Bob's end. This can be studied by considering a unitary coupling 
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between Bob's system B and an ancillary environment system, say C, and then 
tracing C out. Physically, such a quantum operation will emulate environmental 
decoherence. 

We begin by expanding the size of the Hilbert space so that an arbitrary measure- 
ment (or any other quantum operation) can be modeled by coupling to the auxiliary 
subsystem and then discarding it. We assume the ancilla C to initially be in a pure 
state |0), and a unitary interaction U between B and C. Letting primes denote 
the state of the system after U has acted we have H(A, B) = H(A', B'C) as C 
starts out in a product state with AB. We also have I (A : BC) = I (A' : B'C). 
As discarding quantum systems cannot increase the mutual information, we get 
I (A! : B') < I(A' : B'C). 

Now consider the FQSW protocol between A and B in the presence of C. We can 
always view the yield of the FQSW protocol on the system AB to be the same as 
that of performing the protocol between systems A and BC, where C is some ancilla 
(initially in a pure state) with which B interacts coherently through a unitary U . 
Such an operation does not change the cost or yield of the FQSW protocol, as shown, 
but helps us in counting resources. Discarding system C yields 

I {A! : B') < I {A! : B'C) = I (A : BC) = I (A : B). (6.1) 

Now consider a protocol which we call as FQSWD B (FQSW after decoherence), 
where the subscript refers to the decoherence at B. The resource inequality for 
FQSWD B is 

(U S '^ A ' B ' : ^ S ') + h(A' : R')[q -> q] > \l{A' : B')[qq] + (I 5 '~^ : ^ s '). (6.2) 

It is important to explain the terminology used in the above inequality. When we 
have a noisy resource like a mixed state, tjj s ' , or a noisy channel, it is inserted between 
a "()". Thus a mixed state is represented by (vp s '}, and a noisy channel by (Af). A 
channel is a relative resource meaning that the protocol only works 
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provided the input to the channel is the state ijj s ' . On the LHS, U takes the state 
ip s> and distributes it to Alice and Bob. On the RHS, the symbol "I" is an identity 
channel taking the state ijj s ' to Bob alone. The state il) s ' on the left-hand side of the 
inequality is distributed to Alice and Bob, while on the right-hand side, that state 
is given to Bob alone. The primed letters, A', B' etc., indicate that the protocol is 
taking place in the presence of decoherence at Bob's end. 

As in the fully coherent version, Alice is able to transfer her entanglement with 
the reference system R', and is able to distill \l(A' : B') EPR pairs i\qq\) with Bob. 
The net quantum gain for the fully coherent protocol is G = \l(A : B) — \l{A : R) = 
—H(A\B) ebits. This is the difference between the yield obtained and the cost of 
quantum communication incurred. Likewise, the net gain for the protocol suffering 
decoherence at B is G D = \I(A' : B') - \l{A' : R') = -H{A'\B'). Therefore, the net 
advantage of the coherent protocol over the decohered one is given hy D = G — Gd = 
H(A'\B') — H(A\B) ebits. Evidently, this quantum advantage depends on the exact 
nature of the environment and the system's interaction with it via U. Employing the 
original definition of quantum discord due to Zurek [1 1 1] ( Zurek's original definition 
of discord did not consider optimizing over all measurements), D quantifies the loss 
in the yield of a quantum protocol due to environmental decoherence. Our results 
therefore provide a standard way of quantifying, in entropy units, the damage to the 
performance of quantum process and protocols in the presence of any decoherence 
process in any experimental scenario. 

The strength of our result, however, comes from the next step of minimizing D 
over all environmental operations performing measurements. Using the measurement 
model of quantum operations [112], the state pab under measurement of subsystem 
B, changes to p' AB = J2jPjPA\j ® ftj, where {iij} are orthogonal projectors resulting 
from a Neumark extension of the POVM elements [131]. The unconditioned post- 
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measurement states of A and B are 



Pa = ^PjPA\j = Pa, p' B = J2 p ^r 



(6.3) 



j j 



Invoking these relations, we get 



h(A'\b') = J2p,h( Pa \ j ). 



(6.4) 



j 



After minimization over all POVMs, D reduces to T>(A : B) as defined in Eq. (2.29). 
Quantum discord thus quantifies the minimum loss in yield of the FQSW protocol due 
to decoherence. This is our main result, and shows that the performance of all the 
protocols in the quantum information family tree must be judged by the quantum 
discord. The connection between quantum discord and the FQSW protocol provides 
a metric for studying the advantage of coherence in accomplishing any of the children 
protocols that can be derived from the FQSW protocol. For example, we look at 
the noisy versions of quantum teleportation, superdense coding, and entanglement 
distillation. We then apply our general result to quantum computation. 

The connection we have made here is subtle. Although it is known that entangle- 
ment is necessary for the success of the protocol, and that the presence of decoherence 
affects its performance, we have now provided a quantitative result of the amount of 
such a depreciation. We have shown that the amount by which a protocol suffers in 
the presence of decoherence is an inherent property of the quantum states involved. 
It suggests that the best state to be employed in a noisy quantum communication 
protocol is a tradeoff between the entanglement and discord of the state involved. 
The variation of discord and entanglement in quantum states is not monotonic [110], 
and a detailed study of this tradeoff will be presented in future work. In the next 
section, however, we demonstrate the power of our result by applying it to some 
well-known quantum information protocols. 
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6.3 Quantum discord in the children protocols 

The connection between quantum discord and the FQSW protocol provides a metric 
for studying the effect of coherence in accomplishing any of the so called "children 
protocols" that can be derived from the FQSW protocol. In this section, we show 
that by connecting quantum discord with the FQSW protocol, we can interpret 
discord as the advantage of quantum coherence in noisy versions of teleportation, 
super-dense coding, and entanglement distillation. Finally, we reproduce an earlier 
result on the connection of quantum discord and quantum state merging. 

6.3.1 Noisy teleportation 

The noisy teleportation resource inequality can be expressed as 



obtained by combining the mother protocol with teleportation [30]. Here, 1(A) B) = 
—H(A\B) is also known as the coherent information [113]. When Bob undergoes 
decoherence, we get, 



The above can be interpreted as following: The net loss in the number of qubits that 
can be teleported when comparing the coherent teleportation (the one without any 
decoherence), Eq. (6.5), and the one which suffers decoherence, Eq. (6.6), is given by 
I(A)B) - I(A')B') = H(A'\B') - H(A\B). We assume the classical communication 
to be free in this long as we are teleporting unknown quantum states. We 



have H(A\B) = H(A) — I(A : B) = H(A) — I(A : BC) = H(A\BC). As in Sec. (6.2), 



the application of the unitary U, but before discarding the subsystem C, the cost of 



(V AB ) + I (A :B)[c^c]> 1(A) B)[q -> q), 



(6.5) 



(V A ' B ') + I {A' : B')[c ->■ c] > I(A')B')[q q}. 



(6.6) 



teleportation is still given by H(A'\B'C) = H(A\B). From Eq. (6.1), 



H(A'\B') > H(A'\B'C) = H(A\B). 



(6.7) 
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Therefore, we see that the advantage of the coherent protocol over the noisy version in 
teleporting unknown quantum states is equal to the quantum discord of the original 
state. 

6.3.2 Noisy super-dense coding 

Noisy super-dense coding can be derived by combining the mother protocol with 
super dense coding [30] 



showing that one can employ a shared ebit and a single bit of quantum commu- 
nication to communicate 2 bits of classical information. Here, [q — > q] represents 
one qubit of communication between two parties and [qq] represents one shared ebit 
between two parties. Similarly, [c — > c] represents one classical bit of communication 
between the parties. The symbol >z is used to denote exact attainability as com- 
pared to > which is to denote asymptotic attainability. Combining these, the noisy 
super-dense coding protocol can be expressed as, 



When the party B is undergoing decoherence, the noisy superdense coding can be 
expressed as, 



We note that H(A) = H(A'). Thus, due to decoherence, the number of classical bits 
communicated through this protocol gets reduced by the amount I(A : B) — I(A' : 
B'), which is equal to the discord of the original state. 



[qq} + [q^q] h 2[c -> c), 



(6.8) 



(V AB ) + H(A)[q ^q]> I (A : B)[c c]. 



(6.9) 



(V A ' B ') + H{A')[q ^q]> I {A' : B')[c ->• c]. 



(6.10) 
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6.3.3 Entanglement distillation 



The one-way entanglement distillation can be expressed as 



(^ AB ) + I (A :R)[c^c]> I{A)B)[qq\. 



(6.11) 



This inequality can be derived by combining the FQSW protocol Eq. (2.32) and 
recycling the \I{A : R) ebits out of the total hi (A : B) produced for teleportation, 
as shown in [30]. Decoherence at Bob's end B provides 



The net change in entanglement distillation is I(A')B') — I(A)B) = H(A\B) — 
H(A'\B'), which is the negative of the quantum discord of the original state. As is 
well known, classical communication between parties cannot enhance entanglement, 
and we can neglect the overhead of I(A : R) — I{A' : R') classical bits. 

6.3.4 Quantum state merging 

We showed in Chapter 5 that quantum state merging provides an operational in- 
terpretation for quantum discord [4, 139]. It is the markup in the cost of quantum 
communication in the process of quantum state merging, if one discards relevant 
prior information. An intuitive argument for the above interpretation of quantum 
discord can be made through strong subadditivity [114] 



From the point of view of the state merging protocol, the above has a very clear 
interpretation. Having more prior information makes state merging cheaper. In 
other words, throwing away information will make state merging more expensive. 
Thus, if Bob discards system C, it will increase the cost of quantum communication 
needed by Alice in order to merge her state with Bob. 



(* A ' B ') + I {A' : R')[c c] > I(A')B')[qq]. 



(6.12) 



H(A\B,C) < H(A\B). 



(6.13) 
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We can easily derive the results of [ ] starting from the FQSW protocol. We 
start by expressing the quantum state merging as a resource inequality 

<<0 + H(A\B)[q^q) 

+ I (A : B) 9 [c -> c] > (I 5 ^ : ^ s ). (6.14) 

This accomplishes state merging from Alice to Bob at the cost of H(A\B) bits of 
quantum communication. When H(A\B) is negative, Alice and Bob can distill this 
amount of entanglement in the form of Bell pairs. Thus, quantum state merging 
provides an operational interpretation of H(A\B). 

We can derive quantum state merging from the FQSW if the entanglement pro- 
duced at the end of the FQSW protocol can be used to perform teleportation. We 
start by first noting the teleportation protocol [30] 

[qq] + 2[c^c}h[q^q]. (6.15) 

From the FQSW protocol Eq. (2.32), using the entanglement produced at the end for 
quantum communication Eq. (6.15), one gets the quantum state merging primitive 
Eq. (6.14). 

As in Sec. (6.2), the resource inequality for the noisy version of the quantum state 
merging protocol 

(q A ' B ') + H(A'\B')[q^q] 

+ I (A' : B')*[c ->• c] > (I 5 ^ : ^ 5 ). (6.16) 

The cost of quantum communication in this case is S(A'\B'), and the mark up in 
this cost is S(A'\B') — S(A\B), which is equal to the quantum discord of the original 
state. 
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6.4 Quantum discord in quantum computation 

Our approach also allows us to quantify the effect of decoherence on the efficiency of 
quantum computations as a communication protocol. While our work is restricted 
to a particular model of quantum computation (where computation is viewed as 
a communication process), it provides insights into the role of discord in quantum 
computation. In this model, a quantum computation can be described using two reg- 
isters, the memory (M) and the computational register (C), and two programmers, 
the sender and the receiver [115]. At the beginning of the computation, the sender 
encodes the problem to be solved (the message) as a quantum state in the memory 
register. The initial state is operated upon by a sequence of unitaries (the channel), 
and the receiver collects the answer at a later time by looking at the computational 
register C, which contains the output of this communication process when the compu- 
tation is finished. Consequently, the classical capacity of a quantum communication 
channel can be connected to the efficiency of quantum computation. Initially, the 
two registers are uncorrelated. As the computation/communication proceeds, the 
channel, which are the sequence of unitaries that implement the the black box oracle 
queries generates correlations between the two registers. 

Let the memory register of size N be in the state Y2f=i Pj \j)m- The computational 
register is in the state p° c giving the combined initial state Y^j=iPj\j)M <8> Pc- As 



the computation proceeds, the two registers become correlated [J2j=iPj\j) m) <8> 



p c —> Ylij=\Pj\j) m <8> Pc(j)- The mutual information between M and C is given 



by J(M : C) = H(M) - H(M\C) = H(C) - Ej=iPjH (C(j)) , where p M and p c 



are the reduced density matrices of registers M and C respectively and puc is the 
density matrix for the entire system. The mutual information thus calculated is also 
equal to the Holevo bound for the classical capacity of a quantum channel [118]. 
The goal of quantum computation is to maximize this mutual information. Now, 
any quantum operations on C will act to reduce the mutual information between 




117 



Chapter 6. Quantum discord as a resource in quantum protocols 



M and C, i.e. I(M' : C) < I(M : C). Thus, to the degree that I(M : C) is a 
measure of the efficiency of the quantum computation, it is reduced by the amount 
I(M : C) — I(M' : C), which is equal to the discord of state Pmc- 

In a query complexity model of quantum computation, the change of mutual 
information in a single step sets a limit to the scaling of the computation with the 
problem size. If the variation in mutual information I(M : C) in a single step, say 
t, of a computation is \I(M t+1 : C t +i) — I(M t : C t )\ = f(N) bits, then the query 
complexity for the given computation is at least log N/f(N), since logiV bits is the 
maximum mutual information for a system of size N. If the corresponding variation 
in I(M' : C) is \I{M' t+1 : C[ +1 ) - I{M' t : C' t )\ = g(N) bits, then 

\g(N)-f(N)\ < \l{M' t+1 :C' t+1 )-I{M t+1 :C t+1 ) 
-{I{M' t : C' t ) - I(M t : C t )) 
= \V(M t+1 :C t+1 )-V(M t :C t )\ 

< A(M:C) (6.17) 



where A(M : C) is the maximum possible change in quantum discord in a single step 
of the computation. As the query complexity of the computation is now bounded 
by log N/g(N), which in turn is bounded by the consumption of quantum discord 
in the computational process, our result provides a quantitative characterization of 
quantum discord as a resource in quantum computation. 



6.5 Discussions 



The FQSW protocol is the unification of a large class of quantum information theory 
protocols. We developed a vital link between the FQSW protocol and quantum 
discord, providing a unified picture illuminating the role of quantum discord and 
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hence quantum correlations in essentially all quantum communication protocols. 

While all our results are derived for finite-dimensional cases, gaussian quan- 
tum discord [124] has been related to a generalisation of quantum dense coding for 
continuous-variable states, when all the states and operations involved are gaussian. 
The problem was cast as the advantage that can be harnessed by using nonlocal 
quantum interactions. Our results can be extended to continuous variable systems. 

Our work elucidates the role non-classical correlations, as captured by quantum 
discord, play in quantum communication tasks. For an important and large class 
of protocols, quantum discord serves as a metric of how coherently the protocol 
performs. We have been able to quantify the loss quantum communication protocols 
suffer due to decoherence and show that this is aptly captured by quantum discord. 
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Chapter 7 

Summary and Outlook 



7.1 Summary of the main results 

This thesis addresses some fundamental questions in quantum mechanics as seen 
in the light of quantum information theory. We also explore the unique properties 
of quantum systems, like quantum correlations in bipartite quantum systems, that 
serve as a resource and therefore can be harnessed for building future technology and 
engineering applications based upon quantum theory. 

We have studied the relationship between entanglement and chaos for a system of 
isotropically coupled tops in which one of the tops receives a periodic kick around a 
fixed axis. Here the chaos and entanglement arise from the same coupling mechanism 
removing any ambiguity between chaos in the subsystem vs. chaos in the joint-system 
dynamics. The results reported here give further evidence of the fact that chaotic 
systems take quantum initial conditions to pseudo-random quantum states, and that 
the large long-time entanglement average of states undergoing quantum chaotic dy- 
namics is just that of a typical state in the Hilbert space. We see the confirmation 
of this picture in the excellent agreement between the properties of ensembles of 
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quantum states and the numerical results for the eigenvector statistics and long-time 
entanglement average for the completely chaotic system. This approach was also 
found to be highly flexible, applying to subspaces and mixed phase spaces. 

The dynamical generation of random quantum states has implications for the 
dynamical generation of entanglement. It is well known that for large dimensional 
bipartite Hilbert spaces, a random state is highly entangled with almost the max- 
imum entanglement allowed by the dimension [81]. As the large dimensional limit 
is equivalent to the h — > semiclassical limit, and to the degree that the quantum 
analogs of chaotic Hamiltonians generate random states, one expects near maximal 
dynamical generation of entanglement in quantum chaos, to a value that is predicted 
by the statistics at hand. This is not to say that regular dynamics (quantum analogs 
of integrable motion) cannot lead to highly entangled states. Indeed, such behavior 
is seen, and has been previously noted in [62]. Regular dynamics, however, show os- 
cillatory behavior, including in the generation of entanglement. Chaotic dynamics, 
by contrast, lead to quasi-steady state behavior, and typically lead to higher values of 
time-averaged entanglement than regular motion. Taken together, these facts imply 
that the long-time average entanglement in a bipartite system should be a strong 
signature of classical chaos, closely associated with ergodicity in the two dynamical 
descriptions. 

In chapter 4, we found that the rate at which one obtains information about 
an initially unknown quantum state in quantum tomography is correlated with the 
extent of chaos in the system. This is a new quantum signature of classical chaos that 
we have found. In fact, our results can be regarded as signatures of chaos in closed 
quantum systems undergoing unitary evolution. We have been able to quantify 
the information gain using the Fisher information associated with estimating the 
parameters of the unknown quantum state. When the system is fully chaotic, the 
rate of information gain agrees with the predictions of the random matrix theory. 
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In the last two chapters we gave an operational interpretation of quantum dis- 
cord through quantum state merging and established the role of quantum discord 
in quantum communication protocols. Our work elucidates the role non-classical 
correlations, as captured by quantum discord, play in such quantum information 
processing tasks. For an important and large class of protocols, quantum discord 
serves as a metric of how coherently the protocol performs. We have been able to 
quantify the loss quantum communication protocols suffer due to decoherence and 
show that this is aptly captured by quantum discord. In conclusion, this study places 
quantum discord squarely in the midst of quantum informational concepts and opens 
up the way for its manipulation as a resource in quantum information processing. 
We also hope that our work will serve as a stepping stone for a more comprehensive 
and unified understanding of quantum physics, thermodynamics and information 
theory. Our work can be summarized by one slogan: "Quantum discord is the 
advantage of quantum coherence in quantum information theory". 

7.2 Outlook 

One of the indicators of good science is the questions it raises as much as the ones 
it solves. This work can be extended in several interesting directions. I outline 
a few of them. Though quantum systems show no sensitivity to initial conditions 
due to unitarity, they do show a sensitivity to parameters in the Hamiltonian [96]. 
This poses interesting questions for quantum tomography and, more interestingly, 
for quantum simulations. Under what conditions are the system dynamics sensitive 
to perturbations and how does this affect our ability to perform quantum tomogra- 
phy? Under what conditions does the underlying quantum chaos effect our ability 
to accomplish quantum simulations in general? For example, Shepelyansky has done 
extensive work on the issue of many-body quantum chaos in the quantum computer 
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hardware and its effect on the accuracy of quantum computation [23]. Recently 
classical simulations of quantum dynamics have been connected to intergrability and 
chaos [25]. 

Information theoretic characterization of quantum chaos has thrown light on the 
connections between entanglement and chaos [144]. As already mentioned, quan- 
tum chaotic dynamics will drive the system into arbitrary superposition of quantum 
states. For a bipartite system and an initially pure product state of minimum uncer- 
tainty wave packet, this is reflected in the generation of highly entangled states in 
the Hilbert space. What happens when the initial state is a mixed state? I conjec- 
ture that in general, the generation of correlations between subsystems of a chaotic 
system can be explored using quantum discord. Moreover, discord can serve as a 
useful quantity to quantify the amount of quantum information generated when the 
system dynamics are perturbed. 

Our work on chaos has interesting connections to thermalization in closed quan- 
tum systems. We have the evidence that the rate of thermalization in quantum 
systems is closely related to non-integrability, chaos and symmetries respected by 
the system. Could we find a "thermalization witness" based on the measurement 
record and properties of the covariance matrix? Many-body chaos is another avenue 
that can be exlpored using the signatures of chaos we have found. 

Our work on quantum discord raises many interesting questions and opens new 
directions. For example, can a state with zero discord be treated as "classical"? 
Are states with vanishing quantum discord useful for certain tasks which are not 
possible classically? For example, the question whether concordant computations 
can be simulated classically has been investigated by Eastin [100]. A concordant 
computation is one in which after each stage of computation, the resulting quantum 
state is diagonal in a product basis, and hence has zero discord. The entire simulation 
is finding the right product basis, yet Eastin's findings suggest that it might be 
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difficult to simulate such computations efficiently. Therefore, quantum states with 
vanishing quantum discord might have a non-classical character to them. This leads 
to a fundamental question about the border between quantum and classical. Is there 
another way of marking this border? 

Another related direction concerns with the quantification of the amount of non- 
classicality of quantum states. Any measure of non-classicality associates a number 
to a quantum state that quantifies the amount of non-classicality present in the state. 
In this thesis, we have been able to find the operational meaning of that number, 
when the non-classicality was quantified by quantum discord. Is there another way 
of quantifying the quantum character of a system? If the answer is yes, then what 
is the operational significance of such a metric? 



124 



References 



[I] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 777 (1935). 

[2] W. H. Zurek, Annalen der Physik (Leipzig) 9, 855 (2000); H. Ollivier and W. H. 
Zurek, Phys. Rev. Lett. 88, 017901 (2002). 

[3] A. Datta, A. Shaji, and C. M. Caves, Phys. Rev. Lett. 100, 050502 (2008). 

[4] V. Madhok and A. Datta, Phys. Rev. A 83, 032323, (2011). 

[5] Rolf Landauer: "Irreversibility and heat generation in the computing process," 
IBM Journal of Research and Development, vol. 5, pp. 183-191, 1961. 

[6] R. P. Feynman, Rev. Mod. Phys. 20 , 367(1948). 

[7] Steven H. Strogatz, "Nonlinear Dynamics and Chaos" . Westview Press, 1994. 

[8] M. C. Gutzwiller, J. Math. Phys. 11, 1791 (1970). 

[9] F. Haake, Quantum Signatures of Chaos (Spring- Verlag, Berlin, 1991). 

[10] Y. V. Nazarov and Y. M. Blanter, Quantum Transport, Cambridge University 
Press (2009). 

[II] O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev. Lett. 52, 1 (1984). 

[12] T. Bhattacharya, S. Habib, and K. Jacobs, Phys. Rev. A, 67, 042103, (2003). 

[13] R.L. Liboff, Introductory Quantum Mechanics ( Addison- Wesley, Reading, 
1989). 

[14] A. Messiah, Quantum Mechanics (Wiley, New York, 1966) 

[15] Shannon, Claude E. (July/October 1948). "A Mathematical Theory of Commu- 
nication". Bell System Technical Journal 27 (3): 379423. 



125 



References 



[16] W. H. Zurek, Physica Scripta T76, 186 (1998). 

[17] S. Ghose, P. Alsing, Ivan H. Deutsch, T. Bhattacharya, S. Habib, and K. Jacobs, 
Phys. Rev. A, 67, 052102, (2003). 

[18] T. Bhattacharya, S. Habib, and K. Jacobs, Phys. Rev. Lett., 85, 4852, (2000). 

[19] A. Peres, Phys. Rev. A., 30, 1610, (1984). 

[20] R. Schack and C. Caves, Phys. Rev. E, 53, 3257 (1996). 

[21] A. Scott and C. Caves, J. Phys. A: Math. Gen. 36, 9553 (2003). 

[22] Ya.G. Sinai, On the Notion of Entropy of a Dynamical System, Doklady of 
Russian Academy of Sciences, (1959), 124, 768-771. 

[23] K.M. Frahm, R. Fleckinger, and D.L. Shepelyansky, Eur. Phys. J. D 29, 139155 

(2004) . 

[24] J. Emerson, Y. S. Weinstein, M. Saraceno, S. Lloyd and D. G. Cory, Science 
302, 2098 (2003). 

[25] T. Prosen and M. Znidaric, Phys. Rev. E 75, 015202 (2007) 

[26] W. Zurek and J. Paz, Phys. Rev. Lett. 72, 2508 (1994). 

[27] A. Silberfarb, P. Jessen, and Ivan. H. Deutsch Phys. Rev. Lett., 95, 030402, 

(2005) . 

[28] S. Chaudhury, A. Smith, B. E. Anderson, S. Ghose, and P. Jessen. Nature 461, 
768-771 (2009). 

[29] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Infor- 
mation (Cambridge Univ. Press, 2000). 

[30] I. Devetak, A. Harrow and A. Winter, IEEE Trans. Inf. Th. 54, 4587, 2008. 

[31] M. B. Plenio and S. Virmani, Quant. Inf. Comp., 7, 1 (2007); R. Horodecki, P. 
Horodecki, M. Horodecki, and K. Horodecki, Rev. Mod. Phys., 81, 865, (2009). 

[32] T. Cover and J. Thomas, Elements of Information Theory, (Wiley & Sons, New- 
York, 2006). 

[33] A. Abeyesinghe, I. Devetak, P. Hayden, and A. Winter, Proc. R. Soc. A 465, 
2537, (2009). 

[34] M. Mehta, Random Matrices. 2nd edition (Academic Press, San Diego, 1991). 



126 



References 



[35] Cramr, Harald (1946). Mathematical Methods of Statistics. Princeton, NJ: 
Princeton Univ. Press. ISBN 0-691-08004-6. OCLC 185436716. 

[36] J. S. Bell, Rev. Mod. Phys. 38, 447 (1966) 

[37] C.H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. Wootters, 
Phys. Rev. Lett., 70, 1895, (1993). 

[38] C. H. Bennett and S. J. Wiesner Phys. Rev. Lett., 69, 2881, (1992). 

[39] C. H. Bennett, "Quantum information and computing," Physics Today, Vol. 48, 
No. 10, October 1995, pp. 24-30. 

[40] A. K. Ekert, Phys. Rev. Lett., 67, 661, (1991). 

[41] R. Jozsa, N. Linden, Proc. R. Soc. Lond. A, 459, 2011 2032, (2003). 

[42] R. Horodecki, P. Horodecki, M. Horodecki, K. Horodecki, Rev. Mod. Phys. 81, 
865, (2009). 

[43] A. Ekert and R. Jozsa, Rev. Mod. Phys. 68, 733, (1996). 

[44] N. Gisin, G. Ribordy, W. Tittel, H. Zbinden, Rev. Mod. Phys. 74, 145, (2002). 

[45] L. Gurvits, Proceedings of the thirty-fifth annual ACM symposium on Theory of 
computing: 10, 2003. 

[46] M.V. Berry and M. Tabor, Proc. R. Soc. Lond. A 356, 375 (1977). 

[47] M.V. Berry, J. Phys. A 10, 2083 (1977). 

[48] S. W. McDonald and A. N. Kaufman, Phys. Rev. Lett. 42, 1189 (1979). 

[49] H.-J. Stockmann and J. Stein, Phys. Rev. Lett. 64, 2215 (1990). 

[50] L. Henderson and V. Vedral, J. Phys. A: Math. Gen., 34, 6899 (2001). 

[51] For continuous variable systems, the notion of discord, atleast for Gaussian 
states [124], agrees with the lack of a positive-definite Glauber-Sudarshan P 
representation as a nonclassicality criterion, showing that essentially all two-mode 
Gaussian states are a priori nonclassical [138]. 

[52] I. Devetak and A. Winter, IEEE Trans. Inf. Theory, 50, 3183 (2004). 

[53] M. Koashi and A. Winter, Phys. Rev. A, 69, 022309 (2004). 

[54] M. B. Hastings, Nat. Phys., 5, 255 (2009). 



127 



References 



[55] E. J. Heller, Phys. Rev. Lett. 53, 1515 (1984). 

[56] K. Furuya, M. C. Nemes, and G. Q. Pellegrino, Phys. Rev. Lett. 80, 5524 (1998). 

[57] P. Miller and S. Sarkar, Phys. Rev. E. 60, 1542 (1999). 

[58] J. Bandyopadhyay and A. Lakshminarayan, Phys. Rev. R. 69, 16201 (2004). 

[59] Ph. Jacquod, Phys. Rev. Lett. 92 150403 (2004). C. Petitjean and Ph. Jacquod, 
Phys. Rev. Lett. 97 194103 (2006). 

[60] M. Znidaric and T. Prosen, J. Phys. A 36 2463 (2003). 

[61] A. Tanaka, H. Fujisaki, and T. Miyadera, Phys. Rev. E. 66, 45201 (R) (2002). 
H. Fujisaki, T. Miyadera, and A. Tanaka, Phys. Rev. E. 67, 66201 (2003). 

[62] R. Demkowicz-Dobrzanski and M. Kus, Phys. Rev. E 70, 066216 (2004). 

[63] J. Bandyopadhyay and A. Lakshminarayan, Phys. Rev. Lett. 89, 60402 (2002). 

[64] A. Lakshminarayan, Phys. Rev. E. 64, 36207 (2001). 

[65] M. Znidaric, J. Phys. A 40 F105 (2007). 

[66] T. Gorin and T. H. Seligman, J. Opt. B 4 S386 (2002). T. Gorin and T. H. 
Seligman, Phys. Lett. A 309 61 (2003). 

[67] R. F. Abreu and R. O. Vallejos, Phys. Rev. A 73 052327 (2006). 

[68] Y. S. Weinstein and L. Viola, Europhys. Lett. 76 746 (2006). L. Viola, J. Phys. 
A 40, 8109 (2007). W. G. Brown et al, Phys. Rev. E 77, 021106 (2008). 

[69] S. Ghose and B. C. Sanders, Phys. Rev. A 70, 062315 (2004). 

[70] X. Wang, S. Ghose, B. Sanders, and B. Hu, Phys. Rev. E. 70, 16217 (2004). 

[71] I. Percival, J. Phys. B: Atom. Molec. Phys. 6, L229 (1973). 

[72] S. R. de Echaniz et al, unpublished, e-print quant-ph/arXiv:0712.0256. 

[73] R. V. Puri, Mathematical Methods of Quantum Optics (Springer, Berlin, 2001). 

[74] O. Bohigas, S. Tomsovic, and D. Ullmo, Phys. Rep. 223, 43 (1993). 

[75] K. Zyczkowski, Physica E 9 583 (2001). 

[76] T. Gorin, H. Korsch, and B. Mirbach, Chem. Phys. E. 217, 145 (1997). 



128 



References 



[77] W. Wootters, Found. Phys. 20, 1365 (1990). 

[78] F. Dyson, J. Math. Phys 3 140 (1961). 

[79] K. Zyczkowski and M. Kus, J. Phys. A 24 4235 (1994). 

[80] D. N. Page, Phys. Rev. Lett. 71, 1291 (1993). C. M. Caves, unpublished, 
http://info.phys.unm.edu/ caves/reports/measures. pdf . 

[81] P. W. Hayden, D. W. Leung, and A. Winter, Commun. Math. Phys. 265, 95 
(2006). 

[82] O. Giraud, J. Martin, and B. Georgeot, Phys. Rev. A. 76 042333 (2007). 
[83] K. Zyczkowski, J. Phys. A: Math. Gen. 23, 4427 (1990). 

[84] W. Brown, L. Santos, D. Starling, and L. Viola, unpublished, e-print quant- 
ph/arXiv:0707.1331. 

[85] M. Kus, J. Mostowski, and F. Haake, J. Phys. A 21, L1073 (1988). 

[86] S. Zelditch and M. Zworski, Comm. Math 175 1432 (1996). 

[87] E. B. Stechel and E. J. Heller, Ann. Rev. Chem. 35, 563 (1984). 

[88] A. Peres, Phys. Rev. E 53, 4524 (1996). 

[89] J. Emerson et al, Science 302, 2098 (2003). 

[90] R. F. Abreu and R. O. Vallejos, Phys. Rev. A 75 062335 (2007). 

[91] M. V. Berry, Proc. R. Soc. Lond. A 413, 183 (1987). 

[92] S. T. Merkel, C. A. Riofrio, S. T. Flammia, I. H. Deutsch, Phys. Rev. A, 81, 
032126, (2010). 

[93] A. Ben-Israel, T. N. E. Greville, Generalized Inverses: Theory and Applications, 
Springer- Verlag New York, 2003. 

[94] C. A. Riofrio, P. S. Jessen, I. H. Deutsch J. Phys. B: At. Mol. Opt. Phys. 44 
154007 

[95] M. Kus, J. Mostowski and F. Haake, J. Phys. A.: Math Gen 21 L1073. 
[96] A. Peres, Phys. Rev. A., 30, 1610, (1984). 

[97] R. Horodecki, P. Horodecki, M Horodecki and K. Horodecki, Rev. Mod. Phys. 
81, 865 (2007); M. B. Plenio and S. Virmani, Quant. Inf. Comp. 1, 1, (2007) 



129 



References 



[98] E. Knill and R. Laflamme, Phys. Rev. Lett. 81, 5672 (1998). 

[99] A. Datta, S. T. Flammia and C. M. Caves, Phys. Rev. A 72, 042316, (2005); 
A. Datta and G. Vidal, ibid. 75, 042310 (2007). 

[100] B. Dakic, V. Vedral and C. Bruckner, Phys. Rev. Lett. 105, 190502 (2010); A. 
Datta and A. Shaji, Int. Jour. Quant. Inf. 9, 1787, (2011); F. F. Fanchini, L. K. 
Castelano, M. F. Cornelio, M. C. de Oliveira, arxiv: 1106.0289, (2011); B. Eastin, 
arxiv: 1006.4402, (2010). 

[101] L. Roa, J. C. Retamai, M. Alid-Vaccarezza Phys. Rev. Lett. 107, 080401 
(2011), 

[102] B. Li, S. Fei, Z. Wang and H. Fan arxiv:1111.2645, (2011). 

[103] A. Datta, S. Gharibian, Phys. Rev. A 79, 042325 (2009); S. Boixo, L. Aolita, 
D. Cavalcanti, K. Modi, A. Winter, Int. Jour. Quant. Inf. 9, 1643 (2011) 

[104] I. Devetak, A. Harrow and A. Winter, Phys. Rev. Lett. 88, 017901 (2002). 

[105] D. Slepian, J. K. Wolf, IEEE Trans. Inf. Th. 19, 471, (1971); T. Cover and J. 
Thomas, Elements of Information Theory, (Wiley & Sons, New York, 2006). 

[106] J. Oppenheim, arxiv:0805.1065, (2008). 

[107] D. Cavalcanti, L Aolita, S. Boixo, K. Modi, M. Piani, A. Winter, Phys. Rev. 
A 83, 032324, (2011), 

[108] A. Datta, Studies on the role of entanglement in mixed-state quantum compu- 
tation, Ph.D. thesis, University of New Mexico, arxiv:0807.4490, (2008). 

[109] B. Groisman, S. Popescu, and A. Winter, Phys. Rev. A 72, 032317, (2005). 

[110] K. Modi, T. Paterek, W. Son, V. Vedral, and M. Williamson, Phys. Rev. Lett. 
104, 080501 (2010); A. Al-Qasimi and D. F. V. James, Phys. Rev. A 83, 032101 
(2011); D. Girolami, G. Adesso, Phys. Rev. A 84, 052110 (2011). 

[Ill] W. H. Zurek, Ann. der Physik (Leipzig) 9, 855 (2000). 

[112] M. Nielsen and I. Chuang, Quantum Computation and Quantum Information, 
Cambridge University Press (2000), pg 366. 

[113] J. Yard, I. Devetak, and P. Hayden, IEEE Trans. Inf. Th. 54, 3091, (2008). 

[114] M. Horodecki, J. Oppenheim, and A. Winter, Nature 436, 673 (2005). M. 
Horodecki, J. Oppenheim, and A. Winter, Comm. Math. Phys. 268, 107 (2007). 



130 



References 



[115] S. Bose, L. Rallan, and V. Vedral, Phys. Rev. Lett. 85, 5448 (2000). 
[116] P. Shor, SIAM J. Comput. 26, 1484 (1997). (FOCS'94). 

[117] L. Grover, in Proceedings of the 28th ACM Symposium on Theory of Computing 
(ACM, Philadelphia, 1996), p. 212. 

[118] A.S. Holevo, Proble. Peredachi Inf. 9, 177 (1973). 

[119] J. Oppenheim, M. Horodecki, P. Hoodecki, and R. Horodecki, Phys. Rev. Lett. 
89, 180402 (2002). 

[120] W. H. Zurek, Phys. Rev. A 67, 012320 (2003). 

[121] A. Brodutch and D.R. Terno, Phys. Rev. A 81, 062103 (2010). 

[122] S. Luo, Phys. Rev. A, 77, 042303 (2008); L. Mazzola, J. Piilo, and S. Manis- 
calco, Phys. Rev. Lett., 104, 200401 (2010); M. D. Lang and C. M. Caves, ibid., 
105, 150501 (2010). 

[123] J. Maziero, L. C. Celeri, R. M. Serra, and V. Vedral, Phys. Rev. A,80, 044102 
(2009); F. F. Fanchini, T.Werlang, C. A. Brasil, L. G. E. Arruda, and A. O. 
Caldeira, ibid., 81, 052107 (2010); A. Datta, ibid., 80, 052304 (2009); T. Werlang 
and G. Rigolin, ibid. , 81, 044101 (2010); M. Ali, A. R. P. Rau, and G. Alber, 
ibid., 81, 042105 (2010). 

[124] G. Adesso and A. Datta, Phys. Rev. Lett., 105, 030501 (2010). 

[125] W. H. Zurek, Rev. Mod. Phys., 75, 715, (2003). 

[126] M. S. Sarandy, Phys. Rev. A 80, 022108 (2009) 

[127] W. H. Zurek, Phys. Rev. A, 67, 012320 (2003); A. Brodutch and D. R. Terno, 
ibid., 81, 062103 (2010). 

[128] T. Cover and J. Thomas, Elements of Information Theory, (Wiley & Sons, 
New York, 2006). 

[129] M. Horodecki, J. Oppenheim, and A. Winter, Nature, 436, 673 (2005). 

[130] M. Horodecki, J. Oppenheim, and A. Winter, Comm. Math. Phys., 268, 107 
(2007). 

[131] M. A. Neumark, Izv. Akad. Nauk. SSSR, Ser. Mat. 4, 53 277 (1940). 
[132] A. Datta, arxiv:1003.5256 (2010). 



131 



References 



133] I. Devetak, Phys. Rev. A, 71, 062303 (2005). 

134] M. Horodecki et. al, Phys. Rev. A, 71, 062307 (2005). 

135] S. Luo, Phys. Rev. A., 77, 022301, (2008); A. Datta and S. Gharibian, ibid., 
79, 042325, (2009). 

136] M. Berta, arxiv:0912.4495 (2009). 

137] B. M. Terhal, M. Horodecki, D. W. Leung, and D. P. DiVincenzo, J. Math. 
Phys., 43, 4286 (2002). 

138] P. B. Slater, J. Opt. B 2, L19 (2000). 

139] D. Cavalcanti et. al, arxiv: 1008.3205 (2010). 

140] B. Eastin, arxiv: 1006.4402 (2010). 

141] V. Madhok and A. Datta, arxiv:1107.0994 (2011), 6th Conference on Theory of 
Quantum Computation, Communication and Cryptography ( Springer's Lecture 
Notes in Computer Science, 2011). 

142] V. Madhok and A. Datta, submitted to Phys. Rev. A. 

143] A. Datta, V. Madhok, A. Shaji, A. B. Tacla, C. M. Caves, in preparation. 
Preprint available. 

144] C. M. Trail, V. Madhok, and I. H. Deutsch, Phys. Rev. E, 76, 046211, (2008). 
145] V. Madhok, C. A. Riofrio, and I. H. Deutsch, in preparation. 
146] V. Madhok, in preparation. 

147] N. Bohr, Zeitschrift fr Physik 2 (5), 423478, (1920). 
148] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003) 

149] M. Horodecki, J. Oppenheim, and A. Winter, Naure, 436, 673 (2005). 

150] M. Horodecki, J. Oppenheim, and A. Winter, Comm. Math. Phys., 268, 107 
(2007). 

151] V. Cristini, J. Lowengrub and Q. Nie, J. Math. Biol, 46, 191, (2003). 

152] B. Groisman, S. Popescu, and A. Winter, Phys. Rev. A, 72, 032317, (2005); 

153] RW. Spekkens, Phys. Rev. A, 75, 032110, (2007); 

154] N. Harrington and R. W. Spekkens, Found. Phys. 40, 125, (2010); 



132 



